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Quantum-statistical theory of Raman scattering processes 
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Raman scattering from a great number of phonon modes is described from a quantum-statistical 
point of view within the standing-wave model. The master equation for the completely quantum 
case, including laser pump depletion and stochastic coupling of Stokes and anti-Stokes modes, is 
derived and converted to classical equations: either into a generalized Fokker-Planck equation and 
an equation of motion for the characteristic function or into the master equation in Fock repre- 
sentation. These two approaches are developed both in linear and nonlinear regime. A detailed 
analysis of scattering into Stokes and anti-Stokes modes in linear regime, i.e., under parametric 
approximation, is presented. The existence of s-parametrized quasiprobability distributions, in par- 
ticular the Glauber-Sudarshan P-function, is investigated. An analysis of Raman scattering into 
separate Stokes and anti-Stokes modes in nonlinear regime, thus including pump depletion, is given. 
The master equation in Fock representation is solved exactly for the complete density matrix us- 
ing the Laplace transform method. Short-time solutions, steady-state solutions and approximate 
compact form solutions are obtained. Relations between the quasidistribution approach based on 
the Fokker-Planck equation and the density matrix approach based on the master equation in Fock 
representation are presented. The photocount distribution and its factorial moments as well as vari- 
ances and extremal variances of quadratures are calculated in both approaches giving the basis for 
the analysis of the quantum properties of radiation such as sub-Poissonian photon-counting statistics 
and squeezing. A comparison of various statistical moments obtained from numerical calculations 
utilizing our exact solution of the master equation and from the approximate relations for short 
times, as well as obtained under parametric approximation, is presented graphically. 
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I. INTRODUCTION: HISTORICAL DEVELOPMENTS 

Almost simultaneously in 1928 Raman and Krishnan 0, Q and Landsberg 
and Mandel'stamm Q observed a new kind of scattering, now referred to as 
(spontaneous) Raman scattering. For the last 65 years Raman scattering has 
unceasingly been in the forefront of both scientific and experimental investiga- 
tions, particularly after the first observation of stimulated Raman scattering 
by Woodbury and Ng [4[ (see also Ref. [5|). Without exaggeration one can say 
that Raman scattering and spectroscopy constitute a completely autonomous 
discipline. 

The literature on Raman scattering is quite prodigious. The theoretical 
principles and milestone experiments describing the Raman effect are sum- 
marized in a number of excellent monographs and reviews, for instance, by 
Bloembergen @ , Kaiser and Maier Q , Koningstein B , Grasyuk [l(| , Wang 
pll fl2j , Cardona ,Long 14 1 , Hayes and Loudon [15| , Penzkofer et al. [IF 



Kielich [17H19I ] , Shen [20( , D'yakov and Nikitin [2l| , and the most recent reviews 
by Raymer and Walmsley (22[, Pefina [2j|, and Mostowski and Raymer [24j ]. 
We also refer the Reader to the special issue of the Journal of the Optical 
Society of America B (25| which is devoted entirely to Raman scattering. Al- 
though an extensive literature has accumulated dealing with Raman scattering, 
it should be emphasized that the understanding of the fundamental principles 
that govern the process is still incomplete. 

There are several major groups of theories treating the Raman effect in the 
semiclassical and quantum approaches, and theories for standing waves and 
spatially propagating waves. Here, we discuss in detail the quantum theory 
of Raman scattering for several radiation modes only; this implies that the 
theory is the best suited for scattering in a tuned cavity. Nevertheless, some 
predictions from the standing wave model also can be applied for traveling wave 
models [H-[2|. 

Various methods have been applied to the Raman effect in each of the above 
theories. Taking into account the equation of motion as the basis for clas- 
sification, we can distinguish the following approaches, based on the photon 
rate equation, the Schrodinger equation, the Heisenberg equation (Heisenberg- 
Langevin equation), the master equation (generalized Fokkcr-Planck equation), 
and the Maxwell- Heisenberg equation (Maxwell-Bloch equation); we refer to 
Refs. [HI, [U [3(| Hlj]. The above classification is obviously oversimplified. 
Firstly, there are many relations bridging these approaches. For instance, 
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we shall apply the master equation approach from which we shall derive the 
Fokkcr- Planck equation and the photon rate equation. Secondly, there exist 
other alternative methods, which do not fit into our classification. Let us men- 
tion, for example, those developed by Mavroyannis [33 - [33 | and Freedhoff (35j . 
Thirdly, one can classify the Raman effect theories in many other ways (see, 
e.g., Ref. [13). 

We shall be considering the incident laser photons to be scattered by chaotic 
phonons or quantized chaotic vibrations in a crystal. The process leads to 
Stokes and anti-Stokes photons. To the description of Raman scattering, we 
use two trilinear Hamiltonians coupled via an infinite number of phonon modes; 
one Hamiltonian describes Stokes radiation, and the other describes anti-Stokes 
radiation. The problem of coupled Stokes and anti-Stokes modes were been 
studied previously by Bloembergen and Shen [36l - l38j who applied the coupled 
wave theory of nonlinear optics formulated by Armstrong et al. [39|. Later, 
Mishkin and Walls 40] quantized the Stokes and anti-Stokes modes, but dealt 
with the laser mode as a constant amplitude (so-called the parametric ap- 
proximation). In fact, they considered two bilinear Hamiltonians, coupled by 
way of a phonon mode. Stokes scattering was treated as a parametric am- 
plifier, whereas anti-Stokes scattering was treated as a parametric frequency 
converter. A detailed study of quantum statistics of the bilinear Hamiltonians 
proposed by Louisell et al. [4l[, has been extensively carried out e.g., Ref. 
[5Cj an d appli ed to Raman scattering in particular by Mishkin and Walls [4 
Walls |5ll. l52j. Pefina [53-55], Karska and Pefina (5(| and others. Walls [5 
(see also Ref. has extended the bilinear Hamiltonian to a trilinear form 
to describe Raman scattering. The dynamics of Raman processes with trilin- 
ear Hamiltonian has been studied by Szlachetka et al. |58l - l60| . Szlachetka and 
Kielich (6f[, Szlachetka [27j . Trung and Schutte [Hj], Tanzler and Schiitte 63 1, 
Reis and Sharma (64|, Pefina et al. [H, Hfl, Pefina and Kfepelka l67l . 68 j, 
Levenson et al. and others (for general analyses see also Refs. [231 l7Cll. l7lf ) . 

We shall describe Raman scattering from phonons as collective phenomena 
involving the interaction of many molecules. Much attention has also been 
drawn to a microscopic picture of the Raman effect by considering the in- 
teraction with individual molecules. Shen (26l. [72j. in his quantum-statistical 
theory of nonlinear phenomena, proposed the general m + n photon Hamilto- 
nian, describing m emissions and n absorptions, and atomic transitions of an 
ensemble of N /-level atoms. This microscopically correct Hamiltonian con- 
tains Bose operators of a field and Fermi operators for optically active electrons 
and therefore describes a variety of nonlinear phenomena, in particular Raman 
scattering (for two- or three-level atoms). The same general Hamiltonian has 
been used by Walls [52j and McNeil and Walls (73|. Raman scattering from 
a two-level molecular (atomic) system f74T - l84j ] and in a three-level molecular 



system [74 185l - l92l | has been extensively studied by various authors. Walls [74 1 
has shown that a description of Raman scattering from two-level molecules 
with a large cooperation number (coherent molecular coupling) is markedly 
similar to the results for Raman scattering from phonons. This is because for 
coherent molecular coupling sums of the Fermi operators for the individual 
molecules can be replaced by the collective operators approximately satisfying 
boson commutation relations. 

Nonclassical properties of radiation, such as squeezing, sub-Poissonian pho- 
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ton statistics and photon antibunching, remain central topics in quantum op- 
tics. The literature in this area is truly prodigious. The Reader is referred 
to the articles published in Vol. 85 of this series and references therein, for 
instance Refs. [931495I . as well as the reviews b y Kie lich et al. (96|, Leuchs l97l . 
Loudon and Knight [98j , Teich and Saleh [9t| 1 1 OQj ] , Zaheer and Zubairy llOl | 
and the topical issues of the Journ al of the Optical Society of America B [102 [ 
and the Journal of Modern Optics [1031 ] . 

Squeezing properties of Raman scattering have been studied by Pefinova 
and et al. Pefina, Karska and Pefina [56j], Levenson et al. and Pefina and 
Kfepelka [67], [Hj]. Sub-Poissonian photon-counting statistics and/or photon 
anticorrelations (in particular antibunching) have been investigated within vari- 
ous approaches to the standing- wave Raman effec t by Loudon [30j| . Simaan [75|. 
Agarwal and Jh a 18711 . Trung and Schiitte Il04ll Szlachetka and Kielich [6lj . 
Szlachetka et al. [58l - l60j j , Gupta and Mohantv [78l. l79ll . Pefina [53T.I541 ] , Tanzler 
and Schuttejp, Germey et al_Eoi|, Mohanty et aljiq], Krai (ToB| , Gupta, and 
Dash [9(J 1 107) . Ritsch et al. [92| and in papers already mentioned [H, [56|, |65| - 
l68j (see also Ref. [Hj]). We note that photons scattered in the hyper-R aman 
effec t can also exhibit nonclassi cal photo n-counting correlations [§2, fl07l - 
OjJ], H3, HI or squeezing [eil \ll3L [TlH ■ 



We shall analyze, in particular, cross-fluctuations (cross-correlations) in 
quadratures and in photocount statistics between different radiation modes. 
The theory of coherent light scat terin g within the consistent multipole tensor 
formalism, developed by Kielich (ll7l | (see also Refs. [13, US [ill) was success- 
fully applied to disclose a novel cross-fluctuation mechanism. Here, an analysis 
of cross-correlations is presented along the lines of Szlachetka et al. [H, HH (see 
also Ref. [11]), as well as Loudon [30l |. 

We shall be studying sub- or super- Poissonian photon-counting statistics. We 
shall not analyze photon antibunching or bunching. The inclusion, in our Ra- 
man scattering model, of standard (i.e., temporal) photon antibunching would 
pose no problem. Let us mention the differe nce betw een s ub-Poissonian statis- 
tics and anti-bunching pointed out in Refs. [l!8| and |ll9j . which enables us to 
claim that these are distinct phenomena, and definitions should not be con- 
fused. The Raman scattering model is not suitable for in vestigations of spatial 
antibunching as defined by Le Berre- Rousseau et al. |l20| and Bialynicka-Birula 
et al. [l2lj ] in terms of negative angular correlations of photons. 

As stated above, we shall be considering the quantum statistics of Raman 
scattering from phonons. We shall concentrate on a statistical analysis within 
the m aster equation approach to the Raman effect proposed by Shen [26| . 



Walls [1221 ] and McNeil and Walls [73J. This approach has been studie d by 
various author s, e.g ., Simaan [75[, Sch enzle and Brand [3l|, Pefina [H, HUH]; 
Germey et a l. llOEfl , Gupta and Dash Il07ll. Bogolubov et al. [88| , Grygiel [83| , 
Miranowicz [8J] , and Karska and Pefina [56| . Usually a master equation is con- 
verted to a classical differential equation. Here, we shall apply a transformation 
to a Fokker-Planck equation (FPE) for s-parametrized quasidistributions using 
the coherent-state technique and an alternative method of a master equation 
in terms of Fock states (or a rate equation for the conditional photon-number 
probabili ties) . 

Walls [1221 ] was the first to apply the FPE technique to Raman scattering. 
This approach was extensively developed by Pefina and coworkers (Ref. [23| 
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and references therein). Unfortunately, a FPE for the Raman effect has been 
solved exactly under parametric approximation only; i.e., a pump depletion was 
not included. It means that Raman scattering is described as a competing pro- 
cess of parametric amplification (Stokes scattering) and parametric frequency 
conversion of light (anti-Stokes scattering) in a nonlinear crystal. This approx- 
imation seems to be a real shortcoming of the FPE approach. A problem of the 
existence of a solution of the FPE also arises. A diffusion matrix of the FPE 
for the s-parametrized quasidistributions (with s « 1) in many cases is not 
positive or positive semidefinite. Therefore, such a FPE cannot be interpreted 
as an equation of motio n describing the Brownian motion under the influence 
of a suitable force [123] . For this reason the term pseudo- or generalized-FPE 
is used in the literature. It is sometimes argued that equations of this type are 
unphysical. However by doubling the phase space, it is possib le to intr oduce 
a generalized P- representation (the positive P-representation) |l24l . Il25l | . The 
equation of motion for this generalized P-representation is a FPE with a posi- 
tive or positive semidefinite diffusion matrix. The nonpositive definite diffu sion 
matrix plays an essential role in the production of nonclassical fields [126| . 

The second method of an equation of motion in Fock representation has 
been app l ied to various multi-photon Raman processes [30, lZ3, UM, l78H80l l9fA 
Il07l . [Toa IPTl [TT^ - flT6j . The master equation (in terms of Fock states) for 
first-order Stokes scattering can be solved by applying the La plac e transform 
method. Solutions obtained by McNeil and Walls [73|], Simaan [75|], and others 
apply only to the diagonal elements of the density matrix p, which is a serious 
drawback of these formulations [23[ . The photocount statistics (sub-Poissonian 
photon statistics, antibunching, or anticorrelation) can be fully analyzed using 
the diagonal, in Fock represe ntation, matrix elements of p only. However, the 
phase properties of the fields jl27l4l29j | , or squeezing properties [98lfl30t (which 
are sensitive to the phase of the field) require the availability of the non-diagonal 
terms of the density matrix. We shall derive, for Raman scattering including 
depletion of the pump field, an exact solution of the master equation for the 
complete density matrix in Fock representation (n, m\p\n' , m!) with arbitrary 
n, m, n' , m! . 

The classical description of Raman scattering into both the Stokes and anti- 
Stokes fields seems to be well understood [a H3 contrary to quantum de- 
scription, which is hampered by the complexity of the underlying Hamilto- 
nians and hence the complex structure of the equations of motion. One of 
the simplest nontrivial models describing t he co upling of the Stokes and anti- 
Stokes scattering was proposed by Knight [13l|. This Raman-coupled model, 
in which a single atom is coupled to a single-cavity mode by Raman type 
transitions, has attracted some attention and has been generalized to much 
more realisti c ex perimental conditions in t he s ubsequent papers by Phoenix 



and Knight (1321 ] . Schoendorff and Risken 11331 . Agarwal and Puri I134L Il35l | 
as well as Gerry and E berly |l36j . Gerry [l37| . Gerry and Huang ,and 
Gangopadhyay a nd Ray |l39ll . It is quite remarkable that there exists a strict 
operator solution jl34l . Il35| | of the master equation describing the evolution of 
the generalized Knight model, which describes the system of an atom under- 
going Raman transitions between two degenerate levels on interaction with a 
quantized field in a lossy cavity driven by an external field including the effects 
of atomic dephasing collisions. An extension of the propagation theory of Ra- 
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man effect [l4Ct Il4l| to include anti-Stoke s sca ttering has been developed by 
Kilin jl42j and independently by Li et al. jl43| |. As mentioned, we shall ana- 
lyze another model of the Stokes-anti-Stokes coupling wit hin t he framework of 
the temporal theory of Raman effect proposed by Walls |l22| and extensively 
studied by Pefina and coworkers (see Refs. |23j and |68[). 

We shall discuss only temporal variations of fields instead of full tempo- 
ral and spatial analysis. The assumption of monochromatic pump, Stokes, 
and anti-Stokes fields restricts the validity of our theory to a cavity problem. 
However, a temporal evolution in a cavity problem can usually be converted 
to a corresponding steady-state propagation in a dispersionless medium by 
simply replacing the time variable t by —z/c, a "normalized" space variable 
z. This procedure permits us to address nonlinear optical phenomena, in 
particular Raman scattering, in a manner analogous to their classical treat- 
ment [13, HE HE • Formal space-time analogies have also been p ointe d out 
in the differential equations for the propagation of short light pulses |l44| . Ob- 
viously, a full quantum space-time description is considerably more complex 
and resides in solving equations of motion for an infinite number of creation 
(annihilation) operators of the single-mode radiation fields. The total spatially 
dependent field is a sum of the single-mode solutions. 

Here we mention only some spatial propagation theories of Raman scatter- 
ing. For a detailed analysis we refer the Reader to the review by Raymer and 
Walmsley [22| and references therein. The temporal and spatial evolution of 
the radiation fields (laser, Stokes and anti-Stokes fields) in Raman scattering 
was successfully described within the framework of the classical coupled wave 
theory developed by Bloembergen and Shen [361 - 138} (see also Ref. [20]). The 
first quantum theories of Raman scattering including spatial propagation were 
pro pose d independently by von Foerster and Glauber [1451 ] and Akhmanov et 
al. |l46| using the analogy of Raman effect and optical parametric amplific ation 
processes. Another method was propo sed by Emel'yanov and Seminogov [l47l |. 
and Mostowski and Raymer [140L Il4l| usin g the analogy between Raman scat- 
tering and superfluorescent processes |l48| . Spectacular predictions of the lat- 
ter theory have been, in particular, macroscopic pulse-energy fluctuations of 
the emi tted radi ation in a manner reflecting the underlying spontaneous ini- 
tiation jl49l - ll51 |. The ne gativ e-exponential probability distribution (NEPD), 
derived by Raymer et al. jl49l |. describes the macroscopic fluctuations of the 
scattered radiation. Karska and Pefina (5(| pointed out that the NEPD corre- 
sponds to the generating function of the integrated intensity extensively used 
in this paper (see Sect. 5.2). The standing-wave theory of Raman scattering 
properly describes the m acros copic fluctuations in the low-gain and high-gain 
regime (see Refs. [22l l23l . [l43l | and references therein). In the latter limit the 
quantum fluctuations of the generated fields can be thought of as arising from a 
classical noise process, contrary to the low-gain limit, where certain nonclassical 
effects occur. 

Finally we should mention certain crucial experiments revealing some man- 
ifestations of Raman scattering. For more details, see the review article by 
Raymer and Walmsley [22| . Experiments on the detection of fluct uati o ns o f 
Stokes pulse energie s we re carried out by Walmsley and Raymer |l52l Il53| . 
and Fabricius et al. jl54j . The temporal and spatial fluctuations of the Stokes 
beam profile, the spectrum, and delay have been investigated in a number of 
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experiments both for depleted and undepleted pump pulse (for references see 
Ref. [13 )• We mention these experiments because the theory of Raman scat- 
tering for cavity modes, to be presented here, correctly predicts the existence 
of macroscopic quantum fluctuations of the Stokes pulses. 

Generation of Raman solitons in th e he avily depleted pump pulse has recently 
been observed by MacPhe rson et al . jl55j | and Swanson et al. [l56| as predicted 
by Englund and Bowden jl57l Il58| . Cooperative effects in Raman scattering, 
referred to as coo pera tive Raman scattering, which is analogous to two-level 
superfluorescence [l59|, occurs for a las er pu mp not significantly depleted. The 
effect was first observed by Kirin et al. |l 60j j and then re-e xamined under fully 



convincing experimental conditions by Pivtsov et al. [161]. In our analysis we 
clearly distinguish the two cases of the depleted and undepleted pump field, 
and therefore have listed some effects and experiments in which this condition 
for the intensity pump is crucial. 

Hyper- Raman scattering, i.e., the three-p hoton analog of Ram an sc attering, 



was discovered in 1965 by Terhune et al. [1621 ] (see also Ref. 1163 1). This 
eff ect w as predicted theoretically by Neugebauer [1641 ] . Kielich [I65l ll66l] and 
Li [167| prior to its experimental detection. Since the discovery of hyper-Raman 
scattering, numerous papers have appeared reporting theoretical investigations 
and observations of the process in a variety of solids, liquids and vapours. Here, 
we shall not discuss higher- orde r Raman scattering processes. We refer to the 
reviews of Ref s. [l?], H2 H3, [l68] for details and literature. 

This paper is organized as follows. In Section 2, the standing-wave model 
of Raman scattering is constructed and the basic equation of motion (mas- 
ter equation) is derived. In Section 3, we give a short account of multimode 
s-parametrized quasidistributions and s-parametrized characteristic functions. 
In Section 4, we introduce definitions of nonclassical properties of radiation 
such as quadrature ( "usual" and principal) squeezing and sub-Poissonian pho- 
ton statistics. In Section 5, we present the s-parametrized quasidistribution 
formalism of Raman scattering either including (in Sect. 5.1) or neglecting (in 
Sect. 5.2) depletion of the pump laser beam in the process of scattering. In 
Section 6 we develop the density matrix formalism of Raman scattering. We 
derive exact solutions of the master equation in Fock representation in Sect. 
6.1.2. We also give short-time (in Sect. 6.1.1) and long-time (in Sect. 6.1.3) 
solutions of the master equation. In Sect. 6.2 we present approximate solutions 
valid under parametric approximation, i.e., when pump depletion is neglected. 



II. MODEL AND MASTER EQUATION 



Let us analyze Raman scattering starting from a completely quantum Hamil- 
tonian but describing phenomenologically only the net effect, i.e., ignoring the 
details of the scattering mechanism. We describe the interaction of three single- 
mode radiation fields: an incident laser beam at the frequency wl, a Stokes field 
at the frequency ws and an anti-Stokes field at the fr eque ncy wa through an 
infinite phonon system at frequencies wyj, after Walls |l22l ] (see also Refs. [53h - 
(55|), by the effective Hamiltonians: 

Ho = haJLCL^ciL + hujsagas + fiujAa\a A + fr'S^uivj a^- a Vj , (1) 

3 



s 



H s = Uy^XsjaLOsay^ + h.c, 

3 

Ha = fi^^X* A:j a L a\a V j + h.c, (2) 

i 

H T = H a -H s ~ H A , (3) 

where Hs (Ha) is the trilinear interaction Hamiltonian for Stokes (anti-Stokes) 
scattering and H is the unperturbed Hamiltonian. For simplicity, we have 
dropped the zero-point contributions. The annihilation operators for the laser, 
Stokes, anti-Stokes, and phonon fields are denoted by at, as, a a and ayi, 
respectively (we label all Hilbert space operators with caret). The coupling 
coefficient Xsj (Xaj) denotes the strength of the coupling between the Stokes 
(anti-Stokes) mode and the optical phonon at the frequency ojvj- These coef- 
ficients depend on the actual interaction mechanism. In the Hamiltonians ([2} 
we neglect terms describing higher-order Stokes sca tteri ng [57l l74j , as well 
as terms describing hyper- Raman scattering [27l. [65l [l07| . In Section 6.1 in the 
analysis of the Raman effect without parametric approximation we also neglect 
anti-Stokes production. 

In our model we take into account only the radiation modes appropriate for a 
cavity. It should be kept in mind that the several radiation mode description is 
applied to the waves involved in the whole course of the interaction, not only at 
the beginning of the interaction process. This approximation is a shortcoming 
from the experimental point of view, since it is not very suitable for describing 
the most common experimental arrangements used when measuring stimulated 
Raman scattering [20l422l |. 

We apply the rotating wave approximation since in the interaction Hamilto- 
nians © we have omitted terms of the form ayjdgaL + h.c. and ay^a\a,L + h.c.. 
For weak coupling these terms are negligible because they vary rapidly as 
exp[±i(<iVj + \lol — ws,A\)t]i which implies that their average is approximately 
zero for times of evolution much greater than \u)l — ^Js.aI^ 1 , contrary to the 
interaction Hamiltonians Hs and Ha ©j which vary as exp[±i(a;vj — — 
ws,A|)i] giving unity for ujyj ~ \ul — w s,a|- We have also neglected terms of 
the form hy^asaL + h.c. and ciyjCiAUL + h.c. These terms, if included, would 
describe a process in which both the Stokes (anti-Stokes) and laser photons are 
annihilated and created in the scattering act. 

The Hamiltonians describe Raman scattering under t he long wave - 
length approximation, which has several important implications [I45lll6t^ . ll70l |. 
Firstly, we can neglect the intermolecular interactions. Each optical vibrational 
mode of the medium is equivalent to a simple harmonic oscillator. Secondly, 
the optical phonon dispersion is negligible. A typical dispersion curve for op- 
tical phonons, toy(ky), is almost flat for wave vectors ky from the interval 
( — 1/A, 1/A), where A is an optical wavelength. In other words, optical wave 
vectors occupy only a very small volume about the origin of the reciprocal 
lattice. Thirdly, a crystal can be treated as a continuum; thus, from the math- 
ematical point of view, sums over lattice sites can be replaced by integrals over 
a volume of the crystal. This long wave approximation is quite realistic for 
optical processes, in particular Raman scattering. 
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A detailed derivation of the Hami ltonian from first principles has been given 
by von Foerster and Glauber |145| in their quantum propagation theory of 
Raman scattering from phonons. Although we deal with modes in a cavity, 
many aspects of their theory recur in our approach. 

In the case of an unbounded medium the momentum is conserved in the 
interaction, i.e., the sum (difference) of the wave vectors ks (fc^) of the Stokes 
(anti-Stokes) photon and ky of the photon involved in the scattering act is 
exactly equal to the laser light wave vector fcj,, 

k L = k s + k v , k L = kA-k V - (4) 

This means that each laser mode interacts strongly only with phonons having 
a single wave vector (one and only one vibrational mode). This is the require- 
ment of translational invariance. Momentum is no longer strictly conserved for 
interactions in a finite medium, since the introduction of boundaries destroys 
the translational invariance of the medium. The strongest interaction is still for 
those modes, which conserve momentum Q and energy (ojvj ~ I w l ~ ws^l); 
nevertheless, in this case the radiation modes arc coupled to a certain range 
of optical phonons whose wave vectors may not satisfy the adequate condi- 
tions (@| by amounts of the order of the reciprocal of the dimensions of the 



medium 145j . The coupling consta nts A s,-, and \aj contain these momentum 
mismatches via phase integrals [271 . Il22j : 



X S j ~ / exp 



-i ( k L - k s - kvj I -r 



d 3 r, 



v 



V 



Hence, the interaction Hamiltonians © are represented by sums over all optical 
vibrational modes that may scatter into or out of the desired mode. This means 
that the coupling of the radiation fields (in particular Stokes and anti-Stokes) 
through a large number of optical phonons is treated stochastically. 

In the Hamiltonians ([I])-© we have assumed all the radiation fields and 
phonons to be polarized linearly in the same direction. We have not in- 
cluded explicitly the polarizati on states of those p hotons, w hich might affect 
the photon- count ing statistics [l7ll Il72j , squeezing jl73l . Il74| and other prop- 
erties (Ref. [175| and references therein). Obviously, this would require the 
discussion of correlation tensor s, in place of correlation functions, involving the 
photon polarization states [H, Il76l Il77l | . 

The model under discussion is restric ted t o th e approximation of electric- 
dipole transitions. In previous papers jl78l . [l79] , Kielich has proposed and 
extensively developed the formal quantum theory of first-, second-, and higher- 
order processes (in particular Raman scattering) taking into account multipolar 
electric and magnetic quantum transitions. 

A lot of attention has been devoted to a simpler completely boson Hamilto- 
nian applied to the description of the statistical properties of Raman scatter- 
ing by phonons treated as a single monochromatic mode (Refs. [3 [57], HH, l69| 
and [23[ with references therein). It is clear that the use of a large number 



10 

of phonon modes (a phonon bath) in the model Hamiltonians ([2]) provides a 
fuller picture of the scattering processes. In particular, the model describes 
the stochastic coupling of the Stokes and anti-Stokes modes through a phonon 
bath. The assumption of a single phonon mode implies that the Stokes and 
anti-Stokes fields are co upled in a deterministic manner, which seems to be a 
rather serious drawback [l22|. 

As a digression, let us mention that the same phenomenological Hamilto- 
nians @ have been use d in the description of Brillouin scattering (see, for 
example, Refs. jl05l . ll80| and [23[ and references therein). The main difference 
between Brillouin and Raman scattering lies in different kinds of the scatterers 
responsible for these effects: acoustic phonons in the Brillouin effect, and opti- 
cal phonons in the Raman effect. This difference is included in the frequencies, 
the coupling constants Xsj, Xaj and the reservoir spectrum. More important, 
acoustic phonons exhibit much greater dispersion than optical phonons. In our 
approach to Raman scattering we neglect dispersion. This assumption applied 
to Brillouin scattering has considerably less validity. 

We are interested only in the statistical properties of the radiation fields (the 
pump and scattered beams) considered as a system. We therefore remove the 
unnecessary information about the infinite system of optical phonons, treated 
as a reservoir (heat bath). The procedure leading to the master equation is 
widely used in quantum optics. For a general review o f the mas ter equation 
methods and the extensive bibliography see Refs. [23L l5llll8llll82| . We rewrite 
the interaction Hamiltonians H$ and Ha in the interaction picture as: 

4 

H s + H A = nJ2hQk, (6) 

k=l 

where 





= F+ 


= ^2 ^Sj&Yj exp[iu Vj (t - 

3 


to)], 




F 3 


= F? 


= K\j&Vi exp[— iu V j {t 

3 


-to)], 




Qi 


= Qi 


= a L a1; exp[-iQ. s (t - t )], 






Qs 


= Qi 


= a L a\exp[in A (t-t Q )]. 




(7) 



The Qi (Fi) are respectively functions of the system (reservoir) operators only. 
The 'cavity' frequencies il s , ft a are equal to 

tts,A = \ul - ws,a\ ■ (8) 

Since the system and the reservoir variables are mutually independent as it 
follows from 

[a i ,af]=5 ij for i,j = L, S, A, VI, V2, (9) 

we may trace, in standard manner, the complete density matrix over the reser- 
voir leading to the reduced density matrix p(t). Obviously, we cannot obtain 
any reservoir averages from p(t). There are some Raman scattering models (e.g. 
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Refs. (23l . ll05lll45j |). where optical phonons are included in the system, whereas 
other crystal excitation modes, such as acoustical phonons, electric excitations, 
and other species of molecular vibrations, serve as a thermal reservoir. 

The radiation fields are weakly coupled to the thermal reservoir. The anti- 
Stokes mode loses energy to the reservoir. The fluctuations in the reservoir also 
couple back into the system introducing noise from the reservoir. However, we 
apply the Markov approximation, a condition sufficient to ensure that energy 
that goes into the reservoir will not return to the radiation fields. This con- 
clusion follows from the definition of the Markovian system as one that cannot 
develop memory - the future of the system is determined by the present and not 
its past [l8l[ 1 182( j . The importance of this assumption is sometimes stressed in 
the concept of a Schrodinger-Markov (or Heisenberg-Markov) picture, meaning 
the standard pictures under Markov approximation jl8l| . The importance of 
non-Markovian effec ts in Raman scattering has been recently studied by, e.g., 
Sugawara et al. [l83| and Villaeys et al. (9JJ- Obviously, the system operators, 
Qi, obey the same commutation relations under this approximation as they did 
originally. 

To obtain the equation of motion for the reduced density matrix p(t), one 
has to compute the reservoir spectral densities 



w,, = / e lul ' T 



FfFiiT)) dr, (10) 







where (■•■)r is the average over all reservoir operators; Ui takes the values 
±£ls,A- The infinite system of optical phonons is assumed to be densely spaced 
with the number of modes between and u>i + dioi equal to g(uii)diOi, so we 
may replace the sums over the optical vibrational modes by integrals 

oo 

£(...) « / du^-X...). (11) 
i o 

Let us introduce two quantities: AO - the frequency mismatch and - the 
medium "cavity" frequency, defined by 

n = ns+ 2 n \ (12) 

The frequency mismatch Afi, in general, is not equal to zero. It is quite 
realistic for optical phonons to assume that the coupling constants \s,A{wj) 
and the phonon density of g{ujj) are flat in the vicinity of fl, so that we can 
write 

g(n±An) « g(Q), 
\ k (tt±An) « Afc(fi), k = S,A. (13) 
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The reservoir is supposed to be at thermal equilibrium. The phonons are 
unaffected by interaction with the radiation fields. In the classical sense this 
means that the phonons are so quickly damped that they remain in their steady 
state [20I [37l [38| . The mean number of phonons in the reservoir mode at 
thermal equilibrium is defined by the Bose-Einstein distribution 



(n(u V j)) 



( hu} Vj , . 



(14) 



where kg is the Boltzmann constant and T is the temperature of the reser- 
voir. Obviously, as the reservoir temperature approaches absolute zero, the 
mean number of phonons (n(uivj)) tends to zero as well. In Sect. 5.2, we 
analyze Raman scattering in a parametric approximation for a "noisy" reser- 
voir ((h(u!vj)) 7^ 0), whereas we study Raman scattering including the pump 
depletion for "quiet" reservoir ((h(ujvj)) ~ 0). After some algebra one obtains 
from (USD, 

wti = (y + i^s) ((nv) + 1) , 

is 
2 



( 



iAu s ) (n v ), 
w ts = (y -iA^Aj (ny), 
^ 3 + 4 = (y + ^a) ((n v ) + l), 
w ti = (~^y + iAujAS^j ((nv) + 1) , 
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^2£d _ {Aujsa) (nv), 

(2|d + , AwsA ) (( ^) + i) 



^ = K)*. (15) 

All other reservoir spectral densities, in particular the diagonal densities 
(for i = 1, ...,4), vanish. For simplicity we have denoted the mean number of 
phonons at frequency fi by (nv) — (n(Slv))- The gain constant for the Stokes 
mode, 7s, the damping constant for the anti-Stokes model 7_a, and the mutual 
damping constants for both scattered fields, is A, Ias, are 

7fc = 2^(ft)|A fc (ft)| 2 (k = S,A), 

ISA = lis = 27rg(n)X s (n)X* A (Q), (16) 

where g(fi), as earlier, denotes the density of the optical phonon modes (the 
reservoir spectrum) at frequency fl. It is seen that the following simple relation 
between the single and mutual damping constants holds: \isa\ = \lAs\ = 
lAls- The frequency shifts, representing the Lamb shift in the frequency fl « 
ilj, are expressed by the Cauchy principle value, V of the integrals: 

00 „ 

A^-P ^'y do, (k = S,A), 
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A. SA = (Au A s)* = -v i g^)^^^^) ^.. (17) 

The only effect of the Auji is to change slightly the frequency f2, so we neglect 
them. Having calculated the reservoir spectral densities we can write the master 
equation for the reduced density matrix p — p(&L, ds,a A , t) as: 

d 1 

-Q t P = lps{[aLa + Sl pa+as\ + h.c.) 

+ 7j i lA ([a~la A , pa L a~X\ + h.c.) 

+ l;lSAe~ 2 ' lAnAt ([aLa^paLO,^] + [a L a^p, a L a^] ) 

+ ^7Ase 2tAnAt ([a^a A ,pa+a s ] + [a+a A p, a+a s ]) 

-(nv)^Hs ([ata s , [a L a+,p]] + h.c.) 

+ ^jA([a L a\,[a+a A ,p]] + h.c.) 

+ 1SA e- 2lAnAt [aia Sl [aia A ,p]} 

+lAse 2 * AnAt [a L a + A ,[a L a + S7 p]]}. (18) 

The term in 7s represents the amplification of the Stokes mode; the term in 
7a describes the loss of energy from the anti-Stokes mode into the reservoir; 
the "f A s and 7sA~terms represent the stochastic coupling between the Stokes 
and anti-Stokes modes through the reservoir; the remaining terms in (hv)ji 
represent the diffusion of fluctuations of the reservoir into the system modes. 
Eq. (p~8|) describes, moreover, the evolution of the laser beam, i.e., the depletion 
of the laser field, the coupling of the field with the Stokes and anti-Stokes fields, 
as well as the diffusion of the reservoir fluctuations into the laser field. The 
interpretation of the 7s (ja) terms as the amplification (attenuation) of the 
radiation fields is as yet intuitive, but will gain in precision on solution of the 
generalized Fokker-Planck equation. The master equation (fT5|) could have been 
written in more compact form; albeit for purposes of interpretation the above 
form is more convenient. 

The master equation (|18[) . m the particular cas e of p arametric approxima- 
tion, reduces to the equation obtained by Walls |l 22j j and Pefina [5J], and 
reduces to that of McNeil and Walls for Stokes scattering alone but with no 
need for the parametric approximation (73|. Our master equation (|18[) differs 
but slightly in the diffusion terms (ny )^j only from the special case of the 
master equation given by Agarwal jl84j | (see also Ref. [73j]). 

The master equati on may b e solved b y various techniques presented in stan- 
dard textbooks [H, [HI [TH [l85l - [T87j . Here, we apply two methods. We 



convert the master equation to an associated classical equation. On the one 
hand, expressing the quantum equation in s-ordered form one obtains the gen- 
eralized Fokker-Planck equation for the s-parametrized quasi-probability dis- 
tributi on, which can be exactly solved for a class of Ornstein-Uhlenbeck pro- 
cesses [l87l Il88| . On the other hand, one can express the master equation in 
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Fock representation, which can be solved, for instance, by the Laplace trans- 
form method [H [ll III ■ 

In the following Sections we analyze three cases. Firstly, we briefly describe 
coupling of the three quantum radiation fields: the laser, Stokes, and anti- 
Stokes beams. The problem simplifies considerably if one assumes narrow 
quasi-probability distributions. Secondly, we apply the parametric approxi- 
mation, which means that the pump field is treated classically. We include the 
coupling of the Stokes and anti-Stokes field through the phonon bath. Thirdly, 
we separately describe either the laser and Stokes mode or the laser and anti- 
Stokes mode, but include the depletion of the pump laser light. In this case we 
assume the heat bath to be "quiet" . 



III. MULTIMODE s- PARAMETRIZED QUASIDISTRIBUTIONS 



A description of the multimode fields via quasiprobability distributions (qua- 
sidistributions, QPDs) or equiy a lentl y via ch arac teristic functions was firs t 
proposed by Glauber, [HI, EM [l90|, Cahill [l9l| and Klauder et al. [l9^ . 
General ordering theorems have been given by Agarwal and Wolf [1931 ] . The s- 
parametrized single-mode quasidis tribu tions and characteristic functions were 
introduced by Cahill and Glauber [l94| . who extensively studied various ways 
of defining correspondences between the operators and functions. For a re- 
cent review of the multimode s-parametrized functional formalism we refer the 
Reader to Ref. [2^|. Here, we list the basic definitions and properties of the s- 
parametrized multimode quasidistributions and characteristic functions useful 
for our further investigations. 

To solve the master equation l|18p . i.e., the operator equation, we use the 
c-number representations W^({o!fc V) a nd C^({/3k}) of the density operator 
introduced by Cahill and Glauber [194} . These representations not only are 
useful as a calculation tool, but also provide insight into the interrelations 
between classical and quantum mechanics. By virtue of the multimode s- 
parametrized displacement operator 



D (s) ({(3k}) = IP (S) (^) = Il ex P (& a * 



(19) 



where the continuous parameter s belongs to the interval (—1,1), one can 
define the s-parametrized multimode characteristic function as the mean value 
of 



C^m}) = Tr\pD^({P k }) 



(20) 



In our situation involving the three radiation modes laser {k — L) 1 Stokes (S) 
and anti-Stokes (A), the simplified notation in Eqs. flT9l and ([2TJ]) stands for 
({/3fc}) = (Pl, As, Pa)- The Fourier transform of the characteristic function 
C (s) ({M) © readily gives the s-parametrized multimode quasidistribution 
W«(K}), 



W (s) ({"fe})= / C«({/3 fc })exp 



. k 



{otkPl ~ alPk) 



d 2 {/Vtt} ■ (21) 
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For completeness we write the inverse Fourier transform, which enables us to 
determine from namely, 



C (S) ({M)= / W«({a*}) 



I exp 



d 2 {a fc /7r} , (22) 



where integration extends over a k in the following sense: 



d 2 {a k /ir} = J] i" -1 ^"* = JJ d(Rea fc )d(Ima fe ) 

fc=L,S,A k=L,S,A 

or over similarly. The normalization is chosen to satisfy 



W (s \{a k })d 2 {a k /i,} =C^(0) =1. 



(23) 



In the thr ee specia l cases of s = —1,0,1 one recognizes the well-known 
QPDs 0, ll94Hl98j . namely the Q function, the Wigner function, and the 
Glauber-Sudarshan P-function, respectively: 



W({a k }) 
P({a k }) 



({a k }\p\{a k }) = 
w- M W^({a k }), 



W^\{a k }), 



(24) 



with M denoting the number of modes (in our analysis M will be equal to 3, 
2, or 1). One can say that the s-parametrized quasidistribution W( s ) (with 
s from the interval (—1,1)) is a continuous interpolation between the P- and 
Q-functions. The Q-function directly determines antinormally ordered expec- 
tation values, the P-function determines normally ordered averages, and the 
Wigner function can be used directly to calculate the averages of symmetrically 
ordered operators. The following relations hold for any parameter s: 



imp ca*r) 

\ fc / (s) 



Tr 



= / w« ({«*}) II K) mfc («*) w * d2 («*M 



n a/3™* d(~f3*) n « c s ({/3fc}) 



(25) 



{/3fc}=0 



where {f3 k } = 0, in the three-mode case, means that (3l — fls — Pa = 0. The 
generally accepted criterion for the definition of a nonclassical field resides in 
the existence of a positive P-function, i.e., a classical state is one whose P- 
function is n o more sin gular than a (5-function and is nonnegative definite (e.g. 

Refs. (HHEEllliQl). Tllis 

means that the quantum statistical properties of 
the nonclassical field cannot be described completely within the framework of a 
classical probability theory. A detailed discussion of the existence of quasidis- 
tributions ({a k }) for the Raman scattering model under consideration is 
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presented in Section 5.2. The Wigner function always exists as a nonsingular 
function, but may assume negative values, and in this sense is n ot a classical 
probability distribution (nevertheless, as was shown by Stenholm [20l| . exper- 
iments always give a positive Wigner function). The Q- function has the prop- 
erties of a well behaved (bounded, nonnegative and infinitely differentiablc) 
classical probability distribution. 

Let us write down the relation between two s±- and S2-parametrized qua- 
sidistributions: 



W^\{a k },t) = (—?—) /expf 2 —y-\a k -p k \ 

\s1-s2J J \ sx-s 2 ^ 



2 



xW^\{p k },t)d 2 {f3 k /7r}, (26) 

where S2 < s\. It is seen that the quasidistribution W^ 2 -* is given by the 
convolution of W^ 1 ) with the multidimensional Gaussian distribution. The 
analogous relation for characteristic functions (|20l) is simpler and valid for any 
si and S2, 

C W (W,*)= C^)({A},*)expf^iEl^! 2 ) • ( 27 ) 

Even in the case when S\ -parametrized QPDs do not exist, the calculation of 
the expectation values (a +m a n )( Sl ) in si order poses no problem. They can be 
obtained from the corresponding si-parametrized characteristic function C^ Sl ' 
in view of Eq. (|25[) or, equivalently, from a s 2 -parameterized quasidistribution 
yy( S2 ) ; which does exist, by means of the relation 

(Q Hr (a k r) {si) = IU k m k \ (^P a^LJ™ («) 



k 



xW^)({a fc })d 2 WA}, (28) 



where L"(i) is the generalized Laguerre polynomial. Alternatively, to ob- 
tain the Si-ordered moments (a +m a n }/ si } one can use the generalized P- 
representation (positive P- representation) 124Ul26l l202l l203j j . 



IV. PHOTON-COUNTING STATISTICS AND SQUEEZING: 

DEFINITIONS 



To investigate nonclassical phenomena such as sub-Poissonian photon- 
counting statistics or photon antibunching, one needs to know the diagonal 
matrix elements in Fock representation of the density matrix p({ak}) only. We 
start from the probability distribution p(n) of the photon number n in the 
fc-mode field within a given volume V of space at the time t, defined by 



{« fc } 



(29) 
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where rife = \a k \ 2 . The s-parametrized quasidistribution W^ s '({a k }, t) ([2Tj) can 
be readily transformed to the following s-parametrized integrated quasidistri- 
bution (intensity distribution) W^(W,t) by means of the (5-function, 

W«(W,t) = Jw {s \{a k },t)s(^\a k \ 2 -W^d 2 {a k /7r}, (30) 

where the variable W can be interpreted as the integrated intensity. The pho- 
todetection equation gives a connection between the continuous integrated qua- 
sidistribution yy WfW - ^ t) a nd the discrete photon-number distribution first de- 
rived by Mandel [204 1205| . This photodetection equation states that the pho- 
tocount distribution p(n) is the Poisson transform of the integrated quasidis- 
tribution W^'(W,t). A generalized photode tection equation for W^({a k },t) 
or for W (s) (W,i) can be written as [23ll206| 




^^aA d 2 {a k /n} 



d 2 W, (31) 



with L^ /_1 (a;) denoting the generalized Laguerre polynomial. We formally 
identify the photon-number distribution (|29[) with the photocount distribution 
pop . There is some slight difference in their physical interpretation, since the 
former distribution describes the probability of having n photons in the mode 
volume V, whereas the latter distribution describes the probability of detecting 
n photons in the detector volume Vdet, defined by its parameters (sensitive 
area, response time, quantum efficiency, etc.). It can be argued, however, (e.g. 
Refs. [207142091] 'i . that there is perfect physical equivalence between the photon- 



number moments obtained from (|29p and the photocount-number moments 
calculated from (pTTj) . under the assumption of ideal detectors. 

The s-parametrized time-dependent generating function (exp(— AW))( S ), 
defined by the Fourier transform of the s-parametrized quasidistribution 
W^({a k }, t) or characteristic function C (s) ({/? fe }, t): 

(exp(-AW(t))) w = / W {s \{a k },t)e n 



= X -M J C ( S )( {/?fc}; i ) exp^-i^|/3 fc | 2 ^d 2 {/3 fc / 7 r} (32) 



enables us to calculate the photon- number distribution p(n, t) and the s-ordered 
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photon-number moments (n k )^ in a particularly simple manner: 



K«) = ^f" — ( 1 + — ^ ) ( CX P ( -^^rr^ 

' 0) 



(-1)" d n ( s-1 \~ M / ( A 



! d\ n V 2 ; \ r \ 1+ MA 



. (33) 



<W% = (-1)*^ <exp(-AW)) w A=Q . (34) 

Eq. (|33[) takes the simplest form for s = 1. Several parameters are widely 
used in the literature to describe the photon- number statistics, e.g.; the Man- 
del Q parameter, the Fano factor, or the normalized second-order correlation 
function. In our analysis we employ the normalized second-orde r fact orial mo- 
ment of the photon- number operators (or integrated intensity) ( [209t ] , [23| and 
references therein) 



7 ( 2) = (( A ^) 2 )(i) = _ x = (Mn k 1)) _ ^ (35) 

(fik) 2 (rik) 2 (fik) 2 



and its generalization, the normalized p-th order factorial moment of the k- 
th and l-th mode (the normalized two-mode cross-correlation function of p-th 
order) 



( P ) _ (Ki) m 1 _ (h k i(hki - l)...(h M -p+l)) , 
" lA^- 1 - Jf^P !> (■><>' 



where rijy — n^. + fii- The higher-order factorial moments (1361) by comparison 
with the second-order moments (|35[) provide us with further information con- 
cerning the photon-number distributions. In view of the fact that n/y is the 

(2) 

sum of the single-mode photon- number operators, the factorial moment 7^ 
can be written as 

(a) _ (( A "fe') 2 )(i) _ ((An fc ) 2 ) (1) + ((An ; ) 2 ) (1) + 2 (An t An ; ) 

("fc;) 2 (h k ) 2 + (hi) 2 + 2(h k )(hi) 

(2) 

The Mandel Q parameter for the mode fc is equal to 7^ (n^), whereas the Fano 

(2) 

factor F is (7^ (fik) + 1) (the photoemciency r\ of the photodetector is assumed 
to be r] = 1). 

Light with photon-number fluctuations smaller than those of the Poisson 
distribution is called sub-Poissonian (or photon-number squeezed) and is de- 
scribed by a negative value of ^ 2 \ both for ^ 2 ^ in the single-mode case and 
(2) 

for 7^/ in the two-mode case. In Section 6 we analyze the two-mode model of 
the Raman effect that comprises the laser (L) and the Stokes mode (S). We 

show that the sum of photon-number operators in both modes is a constant 

(v) 

of motion, which implies that the factorial moments 7^ are constant as well. 
Henceforth we shall be applying another definition to investigate two-mode 
cross-correlation, referred to as the interbeam degree of second-order coher- 
ence, given by (Ref. [3(| and references therein) 

(2) _ (AhkAhi) _ (fikfn) _ 
9kl ~ (n k )(m) ~ (n k )(n) { ' 
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(our definition deviates from those of Ref. [30( by the extra term —1). 

To investigate squeezing properties of light we introduce the Hermitian 
single- and two-mode operators: 

X k {6) =a k e- ie + a+e ie , (39) 



a jk e~ 



(40) 



where a k i = a k + a;. The operator X k (X k i) for 9 = corresponds to the 
in-phase quadrature component of the fc-th (fc-th and l-th) mode (modes) of 
the field, whereas for 9 = n/2 it corresponds to the out-of-phase component. 
For brevity, we use the notation: X k \ = X k (0), X k2 = X k (ir/2), as well as 
X k n = X k i(0) and X k i 2 = X k i{-K/2). The following commutation rules hold: 



X k i, X k2 — 2i 



(41) 



X kll ,X kl2 \=4i. (42) 

Firstly, we shall discuss in brief the single-mode case. The variances of the 
^-dependent quadrature ([3"9"f are 

((AX k (9)) 2 ) = 2Rc [e-™((Aa k ) 2 )] + ({Aa+ Aa k }), (43) 

which obviously give X k \ and X k2 in special cases. The Heisenberg uncertainty 
relation for quadratures, 

((AX kl ) 2 )((AX k2 ) 2 ) >1, (44) 

lays the basis for the definition of "usual" ("standard") squeezing. The state 
of the field is said to be squeezed if the variance X k i or X k2 becomes smaller 
than unity (in general, smaller than the square root of the right side of the 
uncertainty relation for the quadratures). Equivalently, light whose quantum 
fluctuations in the one quadrature are smaller than those associated with coher- 
ent light (minimizing the uncertainty relation) is called squeezed (in the usual 
meaning). Since, for a given quantum state, the variance (|43[) is still dependent 
on 8, the angle 9 can be chosen in a way to minimize (or maximize) the vari- 
ance. Differentiation with respect to 8 leads to the angles 9 + and 9- for th e 
maximal and minimal variances, respectively, given by the relation j210L l21l| 



where the difference between the angles 9 + and is tt/2. On inserting (|45p 
into (|43|) one obtains the extremal variances 

((AX k± f) = ((AX k (9)) 2 ) 

= ±2|((Aa fc ) 2 )| + ({Aa+,Aa fc }). (46) 
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It is noteworthy that the 6*- dependent variance (|43|) can be expressed in terms 
of the extremal variances 

((AX k (9)) 2 ) = ((Al fc _) 2 )cos 2 (fl-tf_) 

+ ((AX fe+ ) 2 )sin 2 (0-0_), (47) 

which is the equation for Booth's elliptical lemniscate i n polar co ordinates [21 1| . 
The principal squeezing, introduced by Luks et al. [210L 1212} . occurs if the 
minimum variance is less than unity: 

((AX fe _) 2 ) < 1. (48) 

From (|46|) it follows that the principal squeezing requires the fulfillment of the 
condition 

(Aa+Aa fc ) < |((Aa fe ) 2 )|, (49) 

whereas the condition for standard squeezing, in view of (|43[) . is 

min{(Aa+Aa fc )±Re[((Aa fc ) 2 )]} < 0. (50) 

The mathematically elegant formalism of principal squeezing (in particular 
other equivalent conditions for principal squeezing) can be formulated using 
the gener alized Heisenb erg uncertainty relation (the Schrodinger uncertainty 
relation) pTfil . Ell : 

((AX kl ) 2 )((AX k2 ) 2 ) > i({AX fel ,Al fe2 }) 2 + l, (51) 

which includes the Wigner covariance (cross-correlation) of the quadratures 
X k i and X k2 equal to 

({Al fel ,Al fe2 }> =4Im[((Aa fc ) 2 )] • (52) 

For extremal variances ((AX k ±) 2 ) the generalized Heisenberg relation reduces 
to the standard uncertainty relation. 

The generalization of the above definitions for the two-mode case is straight- 
forward. By virtue of the commutator (|4"2")h twice as great as for the single- mode 
case (|4ip . the standard and principal squeezing can be defined, respectively, as 

mm{((AX kll ) 2 ),((AX kl2 ) 2 )} < 2, (53) 

((AX kl _) 2 ) < 2. (54) 

We express the two-mode variances and the Wigner covariances in terms of the 
single-mode moments: 

((AX kl ) 2 ) = ((AX k ) 2 ) + ((Al,) 2 ) + 2{AX k AX l ), (55) 
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({Al wl ,AX H2 }) = ({AX kl ,AX k2 }) + ({ax iu AX 12 }) 

+2(AX kl AX l2 ) + 2(AX k2 AX n ), (56) 

where X k l stands for Xj~i(9) (in particular the quadratures). Relations such 
as ([55"T) and (|56|) for the quadratures hold for the two-mode creation and anni- 
hilation operators a kl . The moments ((X k ) 2 ) and ({AX kl , AX k2 }) are given 
by and (|52p . The remaining cross-correlations have the following form in 
terms of the annihilation and creation operators: 

(AX kl AX n ) = 2Re [(Aa k Aai) + (Aa+Aai)] , 

(AX k2 AX l2 ) = 2Rc [ -{Aa k Aai) + (Aa+Aa,)] , 

(AX kl AX l2 ) = 21m [{Aa k A&i) + (Aa+Aa,)] , 

(AX k2 AX n ) = 21m [(AafeAa/) - (Aa+Aa,)] ■ (57) 

Substituting Eqs. (g3J), (JSU) and ([57]) into dSSJ and dHHJ) one obtains explicit 
depe nden cies of the two-mode quadrature moments on the annihilation opera- 
tors [2ll|. 

Alternatively, the single-mode moments (HBl and (|52")) and the condi- 

tions (|49p and (|50[) for single-mode squeezing can be generalized to a two-mode 
case by simple replacement of a k , X k (9) by a k i, X k i{6), showing complete anal- 
ogy between the single- and two-mode descriptions. In particular, the two-mode 
extremal variances are 

((AX kl± ) 2 ) =±2|((Aa fc; ) 2 )| +{{Aa+,Aa kl }), (58) 

by analogy to (|46|) . 



V. FOKKER-PLANCK EQUATION 

A. Raman scattering including pump depletion 

The master equation (ME) is the quantum equation of motion for operators 
and hence it is possible to solve it directly only for a sm all cla ss of models. As an 
example w e cite the Raman-coupled model of Knight [l3l| and its generaliza- 
tions (Ref. jl35| and references therein). Usually the quantum master equation 
is converted to a classical differential equation. Then, standard methods of 
mathematical analysis can be applied. In this Section we present one of the 
most popular methods: transformation to a generalized Fokker-Planck equa- 
tion (FPE) or equivalently to an equation of motion for character istic functions . 
This method is extensively studied in a number of textbooks [H, llSll 1 1821 . 1 1871 ] 
and consists of performing s-ordering of the field operators in the ME (fT8"]) and 
then applying the quantum-classical number correspondence of coherent-state 
technique. The rules for the transformation of the ME into Fokker-Planck equa- 
tions for th e s-p aram etriz ed quasidistribution W^^o;, a* , ^4) are the following 
(e.g. Refs. [Hi and [216j ): 
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a+A 
Aa+ 



i-4 a 



s± 1 d 

2 da 



(59) 



where ^4 is an arbitrary operator; in particular, A can be the density matrix 
p; A is the classical function associated with the operator A; the parameter 
s takes arbitrary values in the range (—1, 1). If necessary, these rules can be 
applied repeatedly. Similarly, we list the rules of transformation of the master 
equation (|18[) to the equation of mot ion for the s-parametrized characteristic 
function C (s) (/3, (3* , A) (e.g. Ref. dH): 



Ah 
aA 

a+A 
Aa+ 



-^ + S ^)C [S \^P\AI 



(60) 



Applying repeatedly the rule (f59| of transformation to the master equa- 
tion (|18l) and after some lengthy algebra we finally arrive at the general- 
ized Fokker-Planck equation for the s-parametrized quasidistribution = 

W {s) {a L ,as,aA,t): 
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(61) 



which is a general izati on of our former relation for 7,4,75,4 7^ and arbitrary 
parameter s (Ref. |217l |). For brevity, we refer to the generalized Fokker-Planck 
equation [23| simply as the Fokker-Planck equation (FPE). The physical inter- 
pretation of (1611) can be given in the same manner as the interpretation of the 
appropriate terms in the ME (fT5| given in Section 2. The FPE (|6Tj) exhibits 
a highly complicated structure. Nonetheless, the equations of motion for the 
mean values (akCti), {a%ai), (with k,l — L,S,A) can be calculated. In 

particular, we obtain 



-({h L (t)) + {h s (t))+{n A (t))) = 0, 



(62) 



with (hk(t)) = (a*,ak). Eq. (|^|) states that the total mean number of photons 
(in all radiation modes) is a constant of motion. 

The FPE (jM]) contains terms of the form g§-<x,a fe ai>V (s) , gf- ^a k aiyV (s \ 

9§7sl7alr aiW(s) ' wri ere ai,aj,a k ,ai = a L ,a* L ,a s ,a* s ,a A ,a* A . It is seen 
that most components of the drift vector are nonlinear to the third order 
in a, and most components of the diffusion matrix are nonlinear up to the 
second order. It is particularly difficult to solve a differential equation with 
such nonlinear diffusion and drift coefficients. Besides, the FPE (|6"Tj) for 
)Aj( s ) = W^^(aL,as,aA,0 with the parameter s ^ ±1 contains third-order 
derivatives in the terms 1 P—-ip—-^—aiW { ' s \ This could be expected since in 

ooii oaj aafc ' 

En 



23j , for instance in the anharmonic oscillator model (for ref- 
there occur third-order derivatives in the FPEs for the 



many models 
erences see Ref. 
Wigner function (s = 0). 

The corresponding equation of motion for the s-parametrized characteristic 
function C( s >(/3l, fis, 0A,t) can be obtained either from the ME (fT8f by per- 
forming the transformation (JBDJ), or from the FPE (|6"Tj) by means of the Fourier 
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transformation (|22j) with respect to the variables oil, cts, O-A- Finally, we arrive 
at the equation of motion for = C^(Pl, Ps, Pa, t): 
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(63) 



Here, we come upon similar difficulties in the way of obtaining an analyti- 
cal solution of (|63|) as in the FPE case ([Fl]) . inherent in the nonlinearity of 
the coefficients of the terms with first- and second-order derivatives as well as 
the presence of terms with third-order derivatives. Nevertheless, in contradis- 
tinction to W^ s \ the existence of a solution for C^ 1 ^ implies, in view of the 
property (|27l) . the existence of a solution for any other parameter S2- 

In view of the particularly complicated structure of Eqs. (I5T1) and (1551) or 
equivalent equations of motion derived within the completely quantum model 
of scattering into both the Stokes and anti-Stokes modes, it would seem that 
a solution in exact closed form cannot be obtained [23j . It is necessary to 
apply further restrictions or approximations in the model to achieve an ana- 
lytical solution of (j6"Tj) . In Sect. 6.1.2 we present a strict analytical solution 
of the two-mode ME (including pump depletion) in terms of Fock states by 
applying the Laplace transform. In Sect. 5.2 we present solutions of two-mode 
linearized FPEs for W^ s \as, ctA, t) and solutions of equivalent equations of 
motion for C^ s \f3s, (3a, t) in the Raman scattering model under parametric ap- 
proximation. In Appendix A we give the solution of a linearized form of the 
three-mode FPE (IBT1) for W^^aL, as, a a, t) properly describing the evolu- 
tion of the radiation fields valid only on the assumption of small fluctuations 
of the fields around their mean values. There, we restrict our considerations 
to the Q-function (s = — 1) to avoid problems of the existence of the quasidis- 
tribution W^ s \aL, as, aA, t) (particularly important in the case of s close or 
equal to 1) and to simplify the third-order FPE (|6 1 1) to second-order, which 
takes place for s = ±1. Within a similar model of Raman scattering from a 
single phonon mode, Szlachetka et al. (58l - [6l| (see also Ref. and Tanzler 
and Schiitte [63j have solved the equations of motion in the short-time approxi- 
mation up to the second power in time. Within the latter (single phonon mode) 
model Per ina and Kf epelka [67], [68| have obtained approximate solutions using 
the approximation of small fluctuations around a stationary solution. 



Here, to find a solution of the ME (|18l) we apply the parametric approxima- 
tion, so no allowance for pump depletion is included. The trilinear Hamiltonians 
Ha , Hs © can be reduced to bilinear functions as a result of the replacement of 
the annihilation operator fiz,, representing the quantum pump field, by the clas- 
sical complex amplitude of the pump field, cl ■ This approximation effectively 
linearizes our model of Raman scattering. Then, the Fokker-Planck equation 
for the two-mode s-parametrized quasidistribution (as, a a, t) takes the 



B. Raman scattering without pump depletion 
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form 
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(64) 



on applying the rules (|55]l to the ME (fT8"|) with complex classical amplitude ez, 
instead of the annihilation operator cll and transforming the variables otk — > 
exp(— iAflAt)otk and /3k — >■ exp(— iAQAt)/3k with the frequency mismatch Afl 
defined by (|12[) . For brevity, we have incorporated the complex amplitude 
into the coupling constants n$ = 7s|ei| 2 , k a = 7_A|ei| 2 , and ks A = K A g = 
Isa&l- The equation (l64l is a generalization for any p aram eter s (s e (— 1, 1)) 
of the FPE given by Walls for the P-function (s = 1) |l22j and by Pefina for 
the P- and Q-functions (s = ±1) [H, Hij]. If we consider production of the 
Stokes radiation only, neglecting anti-Stokes scattering, then Eq. reduces 
to the s-parametrized FPE obtained by Pef inova et al. |218| | . 

We can interpret the FPE AMI) in the same manner as the ME (fT5|) . The first 
term in (|64l) describes the amplification of the Stokes beam, whereas the second 
term describes the attenuation of the anti-Stokes beam; the third term shows 
the coupling between the Stokes and anti-Stokes fields; the remaining three 
terms account for the noise diffusion from the "noisy" (for nonzero temperature) 
reservoir into the system. Contrary to the former equations of motion (1181) , (|6ip 
and (|60p . we lose all information about the depletion of the laser field. It 
is seen that the FPE (|64j) for any quasidistribution even if related to 

the field ordering s ^ ±1, does not contain third-order derivatives, contrary 
to the FPE (|6ip without parametric approximation. Let us note that (|64p 
describes an Ornstein-Uhlenbeck process |188| | since the components of the 
drift vector are linear and those of the diffusion matrix are constant. Various 
methods have been devel oped for solving the equations of motion for Ornstein- 
Uhlenbeck processes [23l Il87| . For instance, expressing the quasidistribution 
W^(as, a A , t) by its Fourier transform (|22[) with respect to the variables as, 
a A , we obtain the following first-order differential equation for the Fourier 
transform, i.e., for the characteristic function C^(/3s,(3 A ,t): 



+ c.c. 
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+ c.c. 
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8+1 



us \{nv) + 
ka (n v ) + 



2 

1 - s 



K-AS 



({n v ) + ^j0 s p A + c.c^ }c( s )(/3 s ,/3A,i). (65) 



Again, to obtain (|2"2")l . one might use the rules (1501) applied to (fT5)l with 
o-L — > cl. Our further results presented in this Section are mainly based on 
very extensive studies carried out by Per in a l53l . l54j , Pef inova and Per ina [48[ , 
and Karska and Pef ina (5(| (see also Ref. |23|| and references therein). How- 
ever, their solutions of the equations of motion for the Raman effect under 
parametric approximation hold only for quasidistributions W*- 1 - 1 (as, a a, t) or 
W(~ 1%> (as,aA,t) and characteristic functions C^^fis, PA,t) related to nor- 
mal and/or antinormal ordering of the field operators. We generalize their 
results to functions related to s-ordering of the field operators, i.e., to an s- 
parametrized quasidistribution W^^as, a a, t) and s-parametrized character- 
istic function ((3s , Pa, t). Let us use, after Ref. 23], the following simplified 
notation for functions characterizing the quantum noise, i.e., the Wigner co- 
variances and variances as well as the mean values of the annihilation operators 
as and a a- 



B ( k s) (t) 
D k i(t) 

D kl (t) 

C k (t) 
&(*) 



1 



({Aa+(t),Aa k (t)}) 
1 



2 

D ik (t) 



({Aa k (t),A ai (t)}), 



1 



DUt) = --({Aa+(t),A ai (t)}), 

<(Aa fc (t)) 2 ), 
(a k (t)), 



(66) 



where k = S, A and {..., ...} is an anticommutator. Assuming the initial condi- 
tion that the Stokes and anti-Stokes fields are stochastically independent, the 
solution of (|65[) for the s-parametrized characteristic function exists for any 
parameter s and is equal to 

C^((3 s ,p A ,t) = exp{ £ [~Bi s) (t)\(3 k \ 2 

k =S,A 

+ fep|+c.c.)+(M|(t)-cc.)] 

+ [D S A(t)P s P A + D S A(t)PsP A + cc] } , (67) 



where 
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1 - 



B^(t) = lB^-(n v )-^-)\U A (t)\ 2 

^ ) + (nv> + ^)|r 1 (7.| 2 -H/M) 
/),.!(/) = (BW + ^ + ili'jcrsl/ir,!/) 



^SA(t) = Cst7 s (*)V2(t) + CS^(t)Vfir(t), 
Cs(t) = CWj^ + QV^), 
C A (f) = C A U 2 A (t) + C* s Vi(t), 

&(*) = E/ S (t)& + Vs(t)&, 

- M*)£a + V A (t)Cs- (68) 

The solution (|67]l for C( s > ((3s , 0A,t) with any parameter s from (—1,1) is, in 
view of the property (|27p . a straightforward generalization of the solutions 
given by Karska and Pefina [56[ (see also Ref. [23|) for C^^Ps, PA,t) related 
with normal or antinormal field operator ordering. Setting initial values C 'a — 
Cs — 0, which implies that DsA(t) = CU(£) = Cs{t) — for any time t, 
the solution (|6"T|) reduces to that of Pefina [53[ . The time-dependent functions 
Uk(t), Vfe(t) (fc = S, A) appearing in ([55} can be expressed as 



v&(<) 



2 



U s (t) = Q2 + Cy + iAfi) Qi, 

tf A (t) = Q2 - ( Y - iAfi) QI, (69) 



in terms of the auxiliary functions 

exp(PiAi) - exp(P 2 Ai) 



Qi = 



Pi -Pi 



2 



9Qi P 1 cxp(P 1 At)-P 2 exp(P 2 A^) 

g2 - ~dT - ^ ' (70) 



P, 2 = \{^^± [C-^f - 4((Afl)» - i^^AQ)]^}. (71) 

It is seen that for the initial moment of time to the functions V&(to) vanish 
and the Uk(to) are equal to unity, so that the initial Wigner covariances -Dfc,;, 
and -Dfc,; (fc, I — S, A) also vanish as a result of the initial condition of zero 
stochastical correlation between the scattered modes. Let us note that on the 
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assumption of the frequency resonant condition Ail — 0, the functions (|69| 
(fTTj) simplify considerably [53( since Pi = ^(ks — k a) and P 2 = 0. This leads, 
in particular, to the relations 



V s (t) + V A (t) = 0, 
U s (t) + U A (t) = 1 + exp ~ KA Af 



(72) 



To obtain the solution of the FPE we perform the Fourier trans- 

form (HU) of C^(p s ,/3 A ,t), which leads to the s-parametrized quasidistribution 
W^(as, «a, t) in the form 



W( s )(a5,«A,i) 
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{(i (s) )- 



1 



S(*)r+o B 4(ai-G(t))' 
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+P 6 (a 5 - fs(*))(ai ~ &(t)) + cc] } 



(73) 



which generalizes the Karska and Pefina result for antinormal ordering [56| 
going over into s-ordering of the field. The time-dependent functions Ei 
(i = 1,...,6) and have been calculated by Pefinova [53] in her analysis 
of quadratic optical parametric processes. Here, we have the following gener- 
alized s-parametrized functions Ei and occurring in (|73[) : 
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(I/*)) 2 = K^(t)K A s >(t)-2B^(t)B A s '(t)K + (t) 
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+2B i h s \t)C\(t)D SA {t)DsA{t) 
+2B A s \t)C s (t)D* SA (t)D SA (t) - 
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(75) 
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(76) 



The two-mode functions W {s) (a s , a A , t) (d) and C^ifis, f3 A ,t) dSTJ reduce 
to the single-mode functions (ctfc, t) and C' 5 -* (/?&, t) (k — S, A) simply by 
setting either as = Ps — or a a — Pa — 0, implying that the coefficients 
Vs(t), V A {t), Ds A (t), and DsA{t) vanish and, for instance, reduces to 

Contrary to the solution (l67l) for the characteristic function C^ s '(/3s, j3 A ,t), 
the solution (|73|) for the quasidistribution W^(as,ce A ,t) may be absent for 
some s-ordering of the field operators depending on the choice of initial field. 
The condition for the existence of the QPD (IT3"|) . i.e., the existence of the 
Fourier transform ([21} of C (s) ((3 S , Pa, t) ([67|>. is that the function K ( A \t), 
LW(t), ReC A (t) +B A s \t), and 



£W = (^)(t)) 1/2 [HeC s (t) + P^] + (ifW(t)) 



-1/2 



ReC^(t)( J DsA(i)-^(t)) 2 -Si s) (t)| J DsA(t)--Ds^(t)| 2 >0 (77) 



should be positive. 



If any of the four functions K { A sl) (t), L^ Sl > 



and ReCyi(i) + B A Sl \t) (for a particular parameter si) 



>(*), £ (S1) W, 

is not positive definite 
everywhere, the equation of motion (|64j) for the si-parametrized quasidistri- 
bution cannot be interpreted as a FPE describing the Brownian motion, i.e., 
the equation is not a "true" FPE. The quasidistribution W^ Sl ^(as,oeA,t) does 
not exist as a positive well-behaved function; still it does exist as a generalized 
function according to the Klauder-S udars han theorem [l95j |. This property 
is a signature of quantum effects (98l4l00j |. Let us note that it is possible to 
use generalized P-representations (positive P-representations) by d oubli ng the 
phase space, as has been proposed by Drummond and Gardiner [124]. The 
generalized P-representations have been applied success fully to solve mast er 
equations of various nonlinear problems (see, e.g., Ref. [124 Il97t |202| . [203} ). 



This method, if applied to our model, requires us to handle eight real variables 
(not counting time), instead of four. 

For initially coherent Stokes and anti-Stokes fields, i.e., satisfying Cs = Ca = 
Dsa — 0, the rather complicated expressions for (f77|) and ([TBI reduce 
to 



= s«(f)B«(t) 



\D SA (t)f 



(78) 
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{3 \t)B^(t) 



\D S A(t)\ : 



(79) 
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It is seen that, in the case of initially coherent fields, the sufficient condition 
for the existence of W (s) (a s , a A , t) ([731) is on ly tnat the function ([75]) 
shall be positive. One obtains further simplification of the problem under the 
assumption of negligible frequency mismatch (Af2 = 0). The functions B^\(t) 
and DsA(t) (EHJ) now reduce to 

4 S) W ^ -^/-(-^/ + -2-^ + <ftv>/ + ) + ^ > 0, 

K S ~ K A \ K S ~ K A K S - K A J 2 

BiHt) = -^f- ( — — — /- + (nv)f + ) + ^ > o, 
k s - ka \«s - k a J 2 



\D SA (t)\ = V^±/_ ( — + (h v )U + 1 ) , (so) 



with 

, ( K S 

In particular, the Wigner function exists, since 



f = m p(*LJ*te)±\. (81) 



m = \ + \ ^-^l- 1 ^)^ + « A ) + « fl ] > 0, (82) 
4 2 ^5 ka 



< for At > 0. (83) 



contrary to the P-function, which does not exist for t > to, since [5c 

KSKA ( Kg Kj\ \ 1 - 

In general the solution (|T3[) at a given time t exists for parameters s less than 

s < (t) + flW (t) + 1 - ^(4 1 )(t)+4 1) W) 2 -4L(D. (84) 

Assuming that the damping constant 7^4 is equal to the gain constant 7s, or 
equivalently k a — Kg — k, we arrive at 

= - s) 2 + 2(1 - s)(l + 2<n v ))/tAi - a/s a (At) a ], (85) 

which is greater than zero for parameters s less than 
1 ( K Ai + l) 2 



s < 



, ) - + {2(h v )-l(l + 2(n v ) + ^ff + l}^}KAt. (86) 

In particular, £M ([53]) for the P-function and flSJJ) for the Wi gner function 
respectively reduce to 

£ (i) = _(«A*y < forAi>0 



2 

4 ' V v "' ' ' 2 



L (0) = 7 + f (n v ) + l) K At > 0. (87) 
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The condition for s fulfilling > cannot be expressed explicitly in a simple 
form in cases with frequency mismatch AO 0. As another example, let 
us assume that the Stokes and anti-Stokes fields are initially chaotic, which 
mathematically differs from our former example of initially coherent state by 
the presence of nonzero initial coefficients Bj. ^ = (fi c hk) (k = S,A). By 
virtue of the relations (|68|) . the functions B^\t), Ds A (t) for chaotic field are 
the same as for a coherent field with extra terms. Here, the function (|75[) 
is found to be 

= (B { ^\t) + (h chS )\U s (t)\ 2 + (n chA )\V s (t)\ 2 ) 



(fl^t) + (n chA )\U A (t)\ 2 + (n chs )\V A (t)\ 2 

n chs )U s (t)V A (t) + (fi chA )U A (t)V s (t)\ 



\D 



SA\ 



which has a form similar to (f?9"| with the same function B^ s \t) and Ds A (t) 
given by ((5U|) . In the case of equal damping and gain constants we obtain 



lis) 



7AA 



[richs) 



J s({n chA ) + (hchs) - 1) + 7Ai 
1 - 







l + s\ 


(rich a) 


(nv) - 


f 2 J 









+ (rich A) (richS) + 



2 

1 - s 



f 

' 2 
1 - s 



It is seen that the Wigner function always exists, since 
V® = jAt((n chS ) + (rich a) + 1) ((n v ) 
+ (rich A 



1-s 



(89) 



flchS) + g 



2<\ n chS) 



(90) 



whereas the P-function exists only for times shorter than 

At < -((h chA ) + (fichs) + iy 1 {(fichA) 2 [(fi c hs) + ((fiv) + 1) 2 
7 



1/2 



+ (flchA)(flchs)[(flchs) + (flv) 2 + ((flv) + I) 2 ] + (flch s) 2 (flv) 2 } 

+(n chA )((h v ) + 1) + (h chS )(n v ). (91) 



The relation (|89[) is quadratic in s and readily gives an analytic expression for 
the largest parameter s (s < 1) for which the quasidistribution W^^ag, a^, t) 
exists at a given time of the evolution At = t — to- 
la Fig.l we present the function L^ s '(t) for different values of the frequency 
mismatch AO, of the mean number of photons (hy), and of the damping 
(ka) and gain (ks) constants. We assume that the Stokes and anti-Stokes 
modes are initially coherent. Thus, for all discussed cases (Figs. la-Id), the 
condition of a positive definite function L^(t) is sufficient for the existence of 
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FIG. 1: The time and parameter s dependence of the function l/ s '(£), related to the 
existence of the QPD W (s) (a s , a A , t), for: (a) k s = 10 s , k a = 10 10 , \AQ\ = 1 ^ 10 6 
(the surfaces coincide in this range of |Af2|), (hv) = 0; (b) ks — ha = 10 , | Afi| = 1, 
\nv) = -t- 100; (c) k s = M = 10 s , |Afi| = 10 6 , (h v ) = 10; and (d) k s = k a = 10 s , 
|Af2| = 10 6 , (ny) = 0. The Stokes and anti-Stokes fields are initially coherent. The 
dashed lines on the surfaces are depicted for L (s) (t) = 0. 

the corresponding s-parametrized QPD. For clarity, the dashed lines in Fig.l 
are depicted for L^(t) = 0. 

We shall briefly analyze a more general situation, which comprises the above 
cases and others. Let us assume after Refs. [56| and [66{ (for a general analysis 
see Refs. [23[ and [219j |) that the Stokes (k=S) as well anti-Stokes (A) modes are 
initially in squeezed states characterized by complex amplitudes parameters 
rfc, and phases (f>%, superposed with a chaotic field, characterized by the mean 
number of chaotic photons {n c hk)- The initial s-parametrized quasidistribution 
(as, oca, to) is then given by 



-1/2 



k=S,A Kk 

4 s) K-^l 2 -Re(^K-6) 2 )] }, 



(92) 



with 



B i S) = B ( k s \t ) = (coshr) 2 + (n chk ) - 



Ck — Cfc(io) 



-exp(z0 fc )sinh(2r fe ), 



(93) 
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which trivially reduces to the quasidistributions of a pure squeezed state 
{{n c hk) — 0), a coherent state ((h c hk) = r k — 0), or a chaotic state 
(r& = = 0). In Sect. 6.2 we analyze another special case of (|93|) with 
rfc = describing a general superposition of coherent and chaotic fields. 

The Raman effect model under parametric approximation is fully specified 
either by the s- parametrized characteristic function C^ s \Ps,pA,t) (|67|) or the 
s-parametrized quasidistribution W( s )(a;s, a a, t) (|73[) . In particular, by virtue 
of the relations presented in Section 4, one can obtain complete information 
about the photon-counting statistics and squeezing properties of the scattered 
fields from (JB7J) or ([73]). 

One can calculate the photon-counting probability distribution p(n) 
from the quasidistribution (as, ctA,t) or integrated quasidistribution 

W (s) (W,t) O by means of (|3ip . or equivalently from the generating func- 
tion (cxp(-AM / (t))) (s) d32j by virtue of ([33]). We apply the latter method, 
which gives us, after insertion of (a s , a A , t) {73) or C^(/3 S , Pa, t) ([57]) 
into (|32[) . the following time-dependent s-parametrized generating function: 



(exp(-AW0) (s) = A- 2 (4 s) )- 1/2 exp (^j , 



(94) 



where the c[ s ^ (jCj ) are polynomials of the fourth (third) order in A 1 : 

4 



4 S) = E( A ^ 

3=0 



1-S 



C 



0) 



EK 1 



1 - s 



(95) 



Adapting the results of Pefinova and Pefina |48| for the coefficients ctj, b 
(j = 0, 1, ...) occurring in ([95]) . one obtains 



fl0 = [-4M 
+{ 



(CaDsaDas + c.c.Ms 



a i 



«2 



B^DsaDsa 



-(C* S K 



(i) 

A 



C A D 



SA 



C* A D 



SA) 



a 



-\[B { pB^D S A + 2Bfc\D SA 



C* S~1* 7~)* 
S^A U SA 



D S aK-\ZsU 



1, 



-[B^i"^ + 2B^C A D SA + C* s C a Das 



+£s^-]&?& + c.c.} 
+ [5 < — > A], 



(i) 



-p^C£ + D sa Dsa)^s + (B^Dsa + C* s D AS )ZsU 
-(B { PDsa + C* s D* sa ^sCa + c.c.] 
+ [S^A], 



(BU+2BV)\Zs 
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+l(C* s e s + DsaZsU - DsaHsCa + cc) 
a 3 = -\ts\ 2 -\U\*, (96) 



1 



b o - 2 s A ~ s A K + 

~2B { i\c* A D* SA D AS + c.c.)+ 1 -K 2 _ 

- l -[C A (C s D 2 SA + C* s D\ s )+c.c] 

+ [S^A], 
h = 2B'£\K { 2 ) -K + )-2{C s D SA D SA + c.c.) 

+[S < — > A], 
b 2 = 2«+4 1) -X + 

+ [5 < — > A], 

6 3 = 2(4 1) +4 1) ), 

64 = 1, (97) 

where [5 < — > A] stands for the preceding terms albeit with interchanged sub- 
scripts S and A. For brevity, we have omitted the time dependence of Eqs. (|96p 
and (|9"T1). Our formulas (l9~4"|) and (j9"5")) are generalizations of the results given in 
Refs. [23], 48] and 56] for s = — 1 to any s. It is seen that the simplest form 
of (|95|) is for normal ordering of the field operators; hence, here, we use only 

this ordering. Pefinova and Pefina [H| have shown that if the polynomial 
has four single roots Afe = — (j) k , the generating function (exp(— XW(t)))^ 
has the form of the fourfold generating function for Laguerre polynomials 

(exp(-AV^))) (1 ) - ntl + AA.r^exp ("i^^) ' ^ 
The field is described by a superposition of signal components 

Ak = iW - K 1 )- 1 £ ai (-X^) 1 (99) 

1=1 1=0 

and the noise components Afe. With Eq. (l9"5)l available, one obtains [48, 56] the 
following photocount distribution 



4 

p(n,t) = ^2 Il exp 



fc li fc 2 f fc 3» fe 4 I — 1 
fc l+ fc 2 + fc 3 + fc 4 = Tl 



-4, 



1 + A ; 



T - V2 ( - A i \ n nm 

(1 + A,)*«+V2r(fc, + 1/2) fcl U(l + A;)/ 



3G 



and its factorial moments 



4 ,fe 



by applying well-known properties of the generating function of the generalized 
Laguerre polynomials L%(x) to the definition relations (|33[) for p(n,t) with 
s = 1 and to the relations (f34]) for (W fe (t))(i). Much simpler expressions are 
found in the special case when the radiation fields are initially superpositions 
of coherent and chaotic fields. From relations (|93l) with = and (|6"5)) it is 
seen that Cg(t) = CA(t) = DsA{t) — 0. The fourfold generating function (|98[) 
reduces to a twofold generating function in the form of (|98p. where the upper 
limit of the product should be replaced by 2 and the square root in (l+AAfc) -1 / 2 
should be omitted. Then the photocount distribution p(n, t) and its factorial 
moments (W (i))(i) become 

P(M) . <„,,-._ (__|l_ _ _£_) £ (») .iAS-d + a,,-<-» 



(=0 



(102) 



(W n (t)) {1) = £ (") L*-, (-gj , (103) 

where L n (x) — L^(x) and the roots A& and coefficients are [Zsj 

i(4 1J (t) - S«(t))(|a| 2 - IC5I 2 ) - + CC.) 

+^s\ 2 + \U\ 2 )- (104) 

The photon-counting statistics of scattering either into the Stokes or anti- 
Stokes mode can be calculated from formulas (|9"8]l - (|101[) . In the single- mode 
case the moments DsA(t) and DsA(t) vanish, considerably simplifying the 
polynomial Op (|95l) . with coefficients bj, to the form 

£« = J] (A- 2 +2A- 1 B«(t)+4 1 '(t)), (105) 

k=S,A 

with the roots Xi^s.A = — (X 1 )k 

\x,2k = B£ ) {t)T\C k (t)\ (k = S,A). (106) 
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The notation Xis.a and X2S.A, instead of Ai.,2,3,4, emphasizes the dependence 
on the single-mode variables in accordance with the assumption of alternative 
scattering into the Stokes or anti-Stokes mode. Analogously, it is seen that 



A lt!l k = \\^(t)\ 2 T \\C k {t)\-\Cl{t)e k {t) + c.c). 



(107) 



On insertion of (|106l) into twofold functions (|98 |) - (|10ip one immediately obtains 
(exp(-AW fc )) (1) = [(1 + AA lfe )(l + AA 2fe )]- 1 /2 

x exp[~\A lk (l + AAu)- 1 - AA 2fc (l + AA 2 fc) _1 ], (108) 



Pk 



(n) = [(l + A lfc )(l + A 2fc )]-^(i + A 2 - 1 )-™cxp(^- T |i^-- T |^^ 



E 



1 + A" 1 

-T{l + l/2)T{n-l+l/2) V 1 + A" 1 



xL, 



-1/2 



An 



Aifc(l + Aifc) 



L 



-1/2 

n—l 



A 



2k 



A 2 fe(l + A 2 fe) 



(109) 



(wm)d) = mx^Yl 

1=0 
-1/2 



xL 



1 fXik 

T(l + l/2)r(n - I + 1/2) ^A 2fe 

A 2k 
A 2 fc 



A\ k \ 
Aifc 



-/ 



(110) 



To obtain the results of Refs. [23j, [56§ and [218] one should replace \\ k by 
E k — 1 and A 2 fc by F k — 1. In particular, assuming that a scattered (Stokes 
or anti-Stokes) mode is initially in a coherent state (thus C k = 0) the mean 
photon numbers (rife) {k = S,A) are 



(n k (t)) = (W k (t)) {1} = Mt)F + B£\t), 



or explicitly 



(ns(t)) = |e s | 2 exp( Ks Ai) + ((h v ) + l)[exp(/s s At) - 1], 



(111) 



(112) 



(n A {t)) = |a| 2 exp(- Kj4 Ai) + (n v )[l-exp(-K jl At)] > (113) 
whereas the mean-square photon numbers are 

(nt) = (W fc 2 >(i) + (W k ) w 

= Mt)\ 4 + Mt)\ 2 (±B k (t) + l) + 2B 2 k (t)+B k (t). (114) 

Then, the normalized second-order factorial moments (1351) are equal to 
7f (*) = B^^lMt^ + B^it)]' 1 

x{|^(t)r[|4(*)| 2 + 4 1) (*)]~ 1 + 1 }- ( 115 ) 
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Let us proceed to analyze squeezing along the lines presented in Section 4. 
We focus our attention on single- and two-mode squeezed light accordi ng to 
the definition of "usual" squeezing and principal squeezing of Luks et al. [l3Ct 
[H3,[2l3]. Using the definitions §6§ of the functions B { k s) (t), D H (t), D m (t), 
and Ck(t) we readily obtain expressions for the moments of the quadratures 
X k and X ki 



((AX fcl , fe2 ) 2 ) = ±2ReC k (t) + 2B k s \t) 



s, 



(IK 



((AX k± ) 2 } = ±2|C fc (t)| + 2B[ s \t) + s, 



(117) 



({AX kl ,AX k2 }) = 4ImC fc (t), 

(AX kl AX n ) = 2Rc[D kl {t)~D kl (t)], 

(AX k2 AX l2 ) = -2Re [D kl (t) + D kl (t)] , 

(AX kl AX l2 ) = 21m[D kl (t) - D kl (t)] , 

(AX k2 AX n ) = 2lm[D kl {t)+D kl (t)], 



(118) 



where, as usual, k,l = S,A and k ^ I. Thus, the two-mode quadrature vari- 
ances now have the form 



((AX, 
((Al 



SAl) 
SA2) 



= ±2Re [C s (t) + C A (t) + 2D SA (t)} 



B { s\t) + B A s) (t) - 2Re D SA (t) + s 



'A 



(119) 



and the extremal variances are 



({AX S A±f) = ±2\Cs(t) + C A (t) + 2D SA (t)\ 



+2 



B { s s) (t)+B A s \t)-2ReD SA (t) 



(120) 



The single-mode squeezing, defined in standard manner, and the single-mode 
principal squeezing require, respectively, that 



\ReC k (t)\\ 



(k = S,A), 



(121) 



whereas the conditions for the two-mode squeezing are, respectively, 



\Re[C s (t) + C A (t) + 2D SA (t)}\\ (s) (s) 
\C si t) + C A (t) + 2D SA (t)\ \ >B ° V +B - (*)-™*DbaM + : (™) 

Examples of the time evolution of (ns( T )), (n|(r)), (&s(t)), and (a|(r)) are 
given by curves C in Figs. 2, 3, 7, and 8, respectively. We assume that the 
Stokes fields are initially in a coherent state (stimulated Raman scattering) or 
in a vacuum state (spontaneous Raman scattering). The rescaled time r is 
defined by t — s- r = tjs- Anti-Stokes scattering is neglected. The phonon bath 
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is at very low temperature, so we put fiy — 0. In Fig. 9 we present the time 
evolution of the extremal variances ((AXs±(t)) 2 ) for fields initially coherent 
with amplitudes equal to aL — v2, as = and assuming that the heat 

bath is "quiet" (i.e., (ny) = 0). In the model under discussion, the variance 
for the Stokes mode, ((AXs(6, t)) 2 ) (curve C in Fig. 9), is independent of 8, 
i.e., ((AXs+{t)) 2 )= ((Als_(r)) 2 ). Hence, squeezing is not observed if the 
initial Stokes mode is in a coherent state. Even if the Stokes field is initially 
squeezed and 7s > ja (not necessarily 7 a = 0), squeezing will rapidly vanish 
due to strong amplification of this mode, which leads to a strong increase in 
quantum noise 56]. The results of this Section (curves C) are compared with 
the exact solutions (without parametric approximation) derived in Sect. 6.1.2 
(curves A) and the short-time solutions of Sect. 6.1.1. 



VI. MASTER EQUATION IN FOCK REPRESENTATION 

The parametric approximation, applied in the previous Section, introduces 
linearization into our Raman scattering model described by the Hamiltoni- 
ans I!])-©- Here, we shall search for a solution to the nonlinear problem, thus 
including pump depletion. The generalized Fokkcr-Planck equation (I6ip and 
the corresponding equation of motion (|63[) for the characteristic function reveal 
the difficulties to be overcome in the complete analysis of Raman scattering into 
simultaneously both the Stokes and anti-Stokes fields from phonons treated as 
a "noisy" {(fiv) 7^ 0) reservoir. Let us assume that the temperature of the 
medium is low. Under this assumption it is quite reasonable to neglect the 
anti-Stokes scattering (7^ = "fsA = JAS — 0) an d, with regard to Eq. (fT4l) . to 
assume that the reservoir is "quiet" ((ny) = 0). Under these approximations 
the master equation (fTSj) reduces to the simple form [73| : 

dp 1 

■7^ = 2 , pa\a s \ + [a L a^p, a^a s ]) , (123) 

where we have introduced a rescaled time t — > r = jst- Let us denote the 
matrix elements of the reduced density operator p in Fock representation by 

(n L ,n s \P(t)\ n' L ,n' s ) = (n,m\p{r)\ n + i/,m + p) = p n , m {v,p, r), (124) 

where for simplicity we identify til — n, and ns = m; p is the degree of off- 
diagonality for the elements of the matrix p for the Stokes mode, whereas v is 
the degree of off-diagonality for the pump laser mode elements. The master 
equation for the matrix elements (|124[) readily follows from Eq. (I123P and can 
be written as 

d 1 

-z-Pnm{vp.T) = -- [n(m + 1) + (n + v)(m + p + l)]p nm (upT) 
or 2 

+ [(?! + l)(n + fc* + l)77l(77l + /x)] ' Pn+l,m—l 

(vpr). (125) 

The equation (|125p for the diagonal matrix elements p n m(00r) reduces to the 
rate equations of Loudon [3(|, and McNeil and Walls [73]. Simaan (75| (cf. 
Ref. [3fJ) analyzed Raman scattering from a gas of two-level atoms. On the 
assumption that almost all the atoms are in their ground state, the Simaan 



rate equation of Ref. 75[ takes the form of Eq. (|125|) for v = p, = 0. 



40 



A. Raman scattering including pump depletion 



1. Short-time solutions 



Before proceeding to derive an exact solution of (|125|) we shall present 
the short-time solutions calculated with the help of the relation (A(t)) — 
Tr{A[ J o(ro) + /5/(To)AT + /5//(To)(Ar) 2 /2]}, where p"(t ) is found by differentiat- 
ing Eq. (I125P with respect to r. The solutions for the mean (n) and mean-square 
number of photons (h 2 ) in the laser mode up to Ar squared are 

(ti(t)) = (h) - (h)((m) + 1)At 

-[(h 2 )((m) + 1) - (h)((m 2 ) + 3<m) + 2)]^£ (126) 



(n>)) = (h 2 ) - {2(h 2 ) - (n))((m) + l)Ar 

-[2(h 3 ){(m) + 1) - (n 2 )(4(m 2 ) + 13(m) + 9) 



+3(n)((m 2 ) +3(m) +2)] 



(A7 



(127) 

where for brevity we set (n p (ro)) = (h p ) and (m p (ro)) = (to p ) (fc = 1,2,3) 
as well as Ar = t — tq. Then the normalized second-order factorial moment, 
Jj*\t), defined by (j3"5"j) . is equal to 

7i 2) (r) = 7 i 2) + [((^ 2 > 2 - (n 3 )(n))(l + (m» 

+ (n)((n 2 ) - (n))(l + (m) - (to) 2 + (to 2 ))] (n)- 3 (Ar) 2 , (128) 

which reduces to the Simaan result 17511: 



7f(r) = ^ ) + [(h 2 )((h 2 )/(h) + l) 
(n)](n)- 2 (Ar) 2 



(h 3 ) 



(129) 



in the special case in which no scattered photons are excited initially, i.e., 
(to) = (to 2 ) = 0. For the initially coherent Stokes and laser modes the factorial 

moment (|128[) reduces to the simple form ^ = \as\ 2 (At) 2 . 
Our short-time solutions for the Stokes mode are 



(m(r)) = (to) + (n)((m) + 1)At 

(rh 2 )(h) + (m)(3(n) - (h 2 )) + 2(h) - (h 2 ) 



(Ar) 



m 3 ){h) + (m 2 )(7(n) - A(h 2 )) + (m)(U(h) - I2(h 2 ) + (h 6 )) 



8(n 2 



(Ar) 3 



(130) 



(to 2 (r)) = (to 2 ) + (n)(2(TO 2 ) + 3(to) + l)Ar 
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FIG. 2: Time behavior of the mean number of the Stokes photons (m) (solid lines) and 
the laser photons (n) (dashed lines) for the initial fields: (a) \at = v2), \as = 0), 
and (b) \<xt = \/2), \as = V0.2). Numerical results with exact solutions of Sect. 
6.1.2 (curves A); short-time approximation of Sect. 6.1.1 (curves B); parametric 
approximation of Sect. 5.2 (curves C); approximate solutions of Sect. 6.2 (curves D). 
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FIG. 3: Time behavior of the mean-square number of the Stokes photons (m 2 ) (solid 
lines) and the laser photons (h 2 ) (dashed lines) for the same cases as in figure 2. 



2(m 3 )(n) + (m 2 )(9(n) - 4(n 2 )) 

+(rh)(13(h) - 9(h 2 )) + 6(h) - 5(n 2 ) 

+ [2(m 4 )(n) + (m 3 )(21(n) - 14(n 2 )) 

+ (m 2 )(73(n) - 72(h 2 ) +8(h 3 )) 
+{m) (l02(n) - 118(n 2 ) + 21(n 3 )) 

(Ar) 3 



(At) 



+48 n - 60 rT + 13 (ff 



(i 



(131) 



On adding Eqs. (|126j) and (|130j) we note that the sum of the mean number of 
photons in both the Stokes and laser modes is constant (at least up to (Ar) 2 ): 



(n(T)) + (m(r)) = (n) + (rh). 



(132) 



Taking a closer look at Eq. (|125[) . which contains only terms with p nrn 
and Pn+i.m-i, one can draw the more fundamental conclusion that the prop- 
erty (|132p holds for any times, in particular for the steady solutions for r — > oo. 
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Equation (|132[) is a special case of (1621) . Actually, we note in view of the mas- 



ter equation (|123[) that the operator o^(t)o,l(t) + a~g (r)as(r) is a constant of 



, is equal to 
^ 3 Ar 



motion. 

The time evolution of the mean values (n(r)), (to(t)), (n 2 (r)), and (to 2 (t)) 
is shown in Figs. 2 and 3 for initially coherent distributions. Curves B are 
obtained from Eqs. (fHo) . (fm)) . (pO)) and (fl"3T|) . 

The factorial moment 7^ (t), in the case of nonzero (m) 

7i 2) (r) = 7 f - 2((to 2 ) - 2(to) 2 - (m))(n)(m) 

- (m 3 ) (to) 2 (n) - (m 2 ) 2 (m) (n) 

+ (TO 2 )(TO) 2 ((n) 2 +2(n)-(n 2 )) 
-(to 2 ) (to) (2 (n) 2 +2(n) - (fi 2 )) 
-3(TO 2 )(n) 2 + (m) 3 (7(n) 2 + 8(h) - 5(n 2 )) 
+ (TO) 2 (10(fi) 2 +4(n) - 3(n 2 )) 



-3(to) (h) 



whereas in the case when all moments 
can be expressed as 



'(At) 2 



(133) 



(2) 

(for k = 1, 2, ...) are zero, 7^ (r) 



7f(r) 



= 2 7 . 



.2) 



1 



v 6(n d ) - 6(n 2 
,Ar 
3 



2 /(n) 



(n)" 2 ' 



(134) 



To obtain a correct time dependence of the factorial moment (|134[) . it is clearly 
necessary to include in Eqs. (|130[) and (|131[) terms at least up to third order 
in t. An equation similar to (|1 34[) has been obtained by Simaan [75j . In Fig. 
(4a) we compare, in particular, our result for the factorial moments of photon 
number in the Stokes mode calculated with Eq. (|133l) and (|134[) (curves B) with 
that obtained from the exact solution (curves A) of the master equation (|125p 
discussed in Sect. 6.1.2. Our Eq. (IT34| gives much better approximation to the 
exact results than Simaan's formula (33) in Ref. [75j . Analogously in Fig. 5, 
the factorial moments for the laser mode calculated with Eqs. (|128p and (|129[) 
(curves B) are compared, in particular, with the exact solutions (curves A). 
The Eqs. (|133l) and (11 34[) reduce, respectively, to 



7 f(r) = 2M 2 |a s r 2 Ar-(2|as| 4 + 3|a s | 



-3K| 



I"sI 2 |c*l| 2 , 



\a L \ 2 \a s \- 4 



(At) 2 , 



(135) 



7 f (r)=l-^Ar. (136) 

for initially coherent radiation fields. 

The corresponding short-time dependence of the cross-correlation (inter- 
beam) function is 



(n(r)TO(r)) = (h)(rh) + [(h 2 )((m) + 1) - (h)((m 2 ) + 2(to) + l)]Ar 
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FIG. 4: Time behavior of the normalized factorial moments 7g (t) for the Stokes 
mode for the same cases (except for curve D) as in figure 2. 

+ [(fi 3 ){(m) + 1) - (h 2 ) (4(to 2 ) + ll(m) + 7) 



+{n){{m 6 } + 6(to 2 ) + ll(m) + 6)] 



(At) 2 



-[3(n 4 )(2(m) + (to 2 ) + 1) - (n 3 )(3(m) + 4(to 2 ) - 1) 
(n 2 )(68(rh) + 43(to 2 ) + ll(m 3 ) + 36) 

-(n)(66(m) + 47(m 2 ) + 14(to 3 ) + (m 4 ) + 32)1 i^lL. 

J /2 



(137) 



On inserting f|126[) . (|130|) . and f)l 3T|) into the definition (f38|) of the intcrbeam 

degree of second-order coherence, 3£s(t), we obtain the following relation for 
the case when photons are initially present in the Stokes mode: 



(n 2 )((rn) + 1) + (h)((m) 2 - (to) - (to 2 ) - 1) 
-(h) 2 ((m) + l)]((n)(m))' 1 AT 

+ { (ft 3 ) (to) ((to) + 1) + (h 2 )(m)(3(m) 2 - 6(m) - 4(to 2 ) - 5) 

-(n 2 )(n)(3(m) 2 + 5(m) +2) 

+ (h)(m) 2 (to) 3 - 3 (to) 2 + (to) (5 - 3 (to 2 )) 

+4(to 2 ) + (to 3 ) + 4 

-(h) 2 [2(to) 3 - 3(to) 2 - 3(m)((m 2 ) + 2) - 2((to 2 ) + 1) 

(At) 2 



+2(ti) 3 ((to) 2 + 2(to) + l)}(TO)" 2 (n)- 



(138) 



otherwise, for the case (to) = (to 2 ) = (to 3 ) = (to 4 ) = 0, we get 

- (n 2 ) 2 /(h) 



7i 2) + (<r 



.At 



-2(n 2 ) +2(h))(h)- 
-(6(h 4 )(h) 2 + 5(ra 3 )(n 2 )(n) 
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FIG. 5: Same as figure 4, but for the normalized factorial moments 7^ 2 '(t) of the 
laser mode. 

-12(n 3 )(n) 2 - 3(n 2 ) 3 - 22(n 2 ) 2 (n) 
,-4 (At) 2 



+26(n 2 )(n) 2 )(n)- 



12 



(139) 



Assuming that the Stokes and laser modes are initially coherent, Eqs. ()138|) 
and (|139[) reduce respectively to 



7 S(r) = -Ar + (K| 2 - 2|a s | 2 M 2 + \a L \>) \a s \~ A -^.. 



(140) 



7§(r) 



At 



(l-13|a L | 2 -8|a L | 4 ; 



(At) 
12 



(141) 



Equation (|137[) . calculated up to the third order in At, enables us to de- 
termine the relation (I139P correct up to At squared, only. Simaan [75| has 
calculated an expression similar to Eq. (|139l) . Examples of the time evolution 
of g L g (t) for initially coherent fields are presented in Fig. 6. Curves B in 
Figs. 6a and b are calculated with Eqs. (|139[) and (|138l) (including terms up 
to At only). Curve S in Fig. 6a is calculated from the Simaan short-time 
approximate solution (32) of Ref. [75[ . One can compare these results (curves 

(2) 

B and S) with g L g(r) obtained from our numerical calculations utilizing the 
exact solution of the master equation (|123p (curves A) . We note the supremacy 
of our short-time approximation (11411) . 

By analogy to the photon-number moments we calculate, in the short- 
time approximation, the mean and mean square of the annihilation operators, 
( a t ( T )) an d (^fe 2 ( r )) f° r both fields (fc = L, S), as well as the cross-correlation 
functions (aJ(T)dg (t)) and (a,E.(T)a<t (t)) . After some algebra, we arrive at 



(&sM> = (at) + Ka L ){a+)^- 



+((a+ 2 a 2 )(a+) - 2(a+a L ) (a+ 2 a s ) - 3(a+a L )(a+)) 



(At) 



(142) 



At 



(«I(r)) = (a + L ) - (a+)((a+a s ) + 1)— + [{&!){&¥ 



-2(a+ 2 a L )((a+a s ) + 1) + (a+)(3(a+a s ) + 1)] 



(At) 5 



(143) 
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FIG. 6: Same as figure 4, but for the interbeam degree of coherence g L s( T )- Addi- 
tional curve S is calculated with the Simaan short-time approximation (Eq. (32) of 
Ref. [75 



+((SrSi)<ar> - (&la L )(apa s ) - 2(a+a L )(a^)) 



.(At) 



(144) 



(arW) = («r> - (ar)((aja s ) + i)Ar + [<ar><ara£) 

-(a+ 3 a L )((a+a s ) + 1) + (a+ 2 )(3(a+a s ) + 1)]^ 



(145) 
Ar 



+«a+ 3 a 2 }<a+> - ll(a+ 2 a L )(a+) - 6(a+ 2 a L )(a+ 2 a s ) 



+5<a+>(a^a s > + (a+)(a+^) + 4(a+)(a+» 



2 

■s. 
(Ar) 2 
8 



(146) 



<a L (r)a+(r)) = (a L )(a+> + ((a+a 2 L )(a+) - (a L )(a^a s ) ~ 2(a L )(a+)) 



At 



+«afa!)(a+) - ll(a+a|)(a+) - 6(a+a| ) (a+ 2 a 5 ) 
+9(a L )(a+ 2 a s ) + (aL)(a+ 3 a|) + 12(a L ) <a+»^. 



(147) 



For brevity, here, we shall restrict our considerations to initially coherent 
states for the Stokes mode denoted as as — \as\exp(icj)s) and for the laser 
mode Qi = |q!l| exp(i^i). In Figs. 7 and 8 we demonstrate the evolution 
of our short-time approximations for (ag(r)) (solid line B in Fig. 7), (cl l (t)) 
(dashed line B in Fig. 7), (o^(t)) (solid line B in Fig. 8), and (o^(t)) (dashed 
line B in Fig. 8) for initially coherent radiation modes. 

From the general relations (|130[) and (|126[) . under the condition of initial 
coherent Stokes and laser fields, we get 

(m(r)) = M 2 + M 2 (M 2 + l)Ar 



4G 




V Of) | N| 
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X 

FIG. 7: Time dependence of the expectation values of the field amplitudes (ag(r)) = 
(as(r)} (solid lines) and (&l(t)) = (ql(t)) (dashed lines) for fields initially coherent 
\oil = \/2} and \cts = V0.2). Curves A, B, C are calculated within the formalisms of 
Sects. 6.1.2, 6.1.1 and 5.2, respectively. 

+ \a L \ 2 [(\a L \ 2 + l)(|a s | 2 + 1) - (|«5| 4 + 4|a s | 2 + 2)]^! (148) 



(n(r)) = |a L | 2 + \a s \ 2 - (m(r)). (149) 

Inserting ([1481) as wel1 as dj) and (fT44"l) with (a^ p a q k ) = \a k \ p+q exp{i(p-q)<j) k } 
(k = S, A), into (j43|) we obtain the ^-dependent variance for the Stokes mode 

((AX s (e)) 2 ) = l + 2\a L \ 2 AT+\a L \ 2 {\a L \ 2 

-[1 + co S 2 (6 - q> S )]\a S \ 2 - 1}(At) 2 . (150) 

The minimal variance ((AA"s_) 2 ), which follows from (|46|) . or directly 
from (11501) . is equal to 

((AXs-) 2 ) = 1 + 2|a L | 2 Ar + \a L \ 2 (\a L \ 2 - 2|a s | 2 - l)(Ar) 2 . (151) 
Analogously, for the laser field we obtain the following ^-dependent variance: 

((AX L (6)) 2 ) = 1 + l[cos(2fl - 20 L ) + l]|a L | 2 |a s | 2 (Ar) 2 , (152) 

on insertion of Eqs. (fT4"9")l . (fT4"51l . and (H"4"5)) into (g21). With regard to the 
relation (|4"6")l , the minimal variance for the field, ((AA"l_) 2 ), is constant up to 
the second order in time: 

((AX L _) 2 ) = 1. (153) 

One can readily deduce the maximal variances ((AXr,,s±) 2 ) from (|150[) 
and (|152[) or from (|46[) . The time evolution of the single- mode extremal vari- 
ances obtained from (I150|) - (|153[) is presented in Figs. 9 and 10 (for (f>L = 0): 
((AaV(t)) 2 ) = ((AA^r)) 2 ) (solid line B in Fig. 9), ((AX s+ (r)) 2 ) = 
((Ai si (r)) 2 ) (dashed line B in Fig. 9), ((AX L _(r)) 2 ) = ((AX l1 (t)) 2 ) (solid 
line B in Fig. 10), ((AX l+ (t)) 2 } = ((AX L2 (t)) 2 ) (dashed line B in Fig. 10). 
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FIG. 8: Time dependence of (a|(r)} = (a|(r)) (solid lines) and (a|(r)} = (a|(r)) 
(dashed lines) for the same cases as in figure 7. 

The covariances for quadratures in the Stokes and laser mode, according 
to (|52| . are, respectively, 

({AX S1 ,AX S2 }) = M 2 |a s | 2 sin(20 s )(Ar) 2 , (154) 

({AX L1 ,AX L2 }) = -|a L | 2 |a s | 2 sin(20 L )(Ar) 2 . (155) 

The generalized Heisenberg uncertainty relation (|5"Tj) with the covariances (|154p 
and (|155p inserted, takes the following form for the Stokes mode in our short- 
time approximation: 

4|a L | 2 AT + M 2 (6|a L | 2 - 3|a s | 2 - 2)(Ar) 2 > 0, (156) 

and for the laser mode 

2|a L | 2 |a s | 2 (Ar) 2 > 0. (157) 

To obtain the two-mode variances and covariances of the quadratures one 
has to calculate, apart from the single- mode functions (I150p . (|152p , (|154[) . and 
(|155[) . the cross-correlations (|57l) . which are obtained in the following form: 

(AX L1 AX S i) = ~\a L \ |as|{2cos0L cos0 s Ar + [cos(0l - <j>s) 

x(4n L - 6n s - 11) + cos(0 L + S )(4|« L | 2 - 2|a s | 2 - 3)^-}, (158) 

(AX L2 AX S2 ) = -\a L \ |a s |{2sin0 L smcf) S AT + [cos(</> L - <f> s ) 

x(4n L - 6n s - 11) - cos(</) L + s )(4|a L | 2 - 2|a s | 2 - 3)^^}, (159) 

(Al L1 Ai S2 ) = | ai,| | a s | {2 cos <^ sin S At + [sin(<^ s - <f> L ) 

x(4n L - 6n s - 11) + sin(0 s + M{A\a L \ 2 - 2\a s \ 2 ~ 3)^^}, (160) 
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FIG. 9: Time dependence of the extremal variances ((AXs-) 2 } (solid lines) and 
((AXs+) 2 ) (dashed lines) for the same cases as in figure 7. 




FIG. 10: Time dependence of the extremal variances {(AXl-) 2 } (solid lines) and 
((AXl+) 2 } (dashed lines) for the same cases as in figure 7. 

(AX L2 AX S i) = |a,s I {2 sin 0l cos <f> s At + [sin(<p L - <f>s) 

x(4n L - 6n s - 11) + sin(<^ + 0s)(4|a L | 2 - 2|a s | 2 - 3)^^}. (161) 
Thus, the two-mode Wigner covariance ([56]) of the quadratures X^si and X^S2 



({AX LS i,AX LS2 }) = 4|a L ||as|sin(0 L + 0s)At 

1 



1 



4|«l| |a s | cos-(0s +^i)sin-(0s - </>l) 

+\a L \\a s \ sin(0 L + <fe)(4|a L | 2 - 2|a s | 2 - 3)] (At) 2 , 

(162) 



The two-mode variances (|55|) of Xlsi and XlS2 are 
((AX LSh2 ) 2 ) = 2 + 2{|a i | 2 -|a L ||as|[cos(0 L -0s)±cos(0 L +0s)]}Ar 
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-{2M 2 (M 2 



a s \ 2 - 1) ± |a L | 2 |a s | 2 (cos20 L - cos20 s ) 



2i |2 

I as 

|a s |M[cos(> L - 4>s){^\a L \ 2 - 6|a s | 2 - 11) 

(At) 2 



± cos(0 L + s )(4|aL| 2 -2|a s | 2 -3)]} 
whereas the extremal variances are 



((AX LS ±f) = 2 + 2[|a L | 2 -|a L ||a s |cos(0 L -0 s )] At 
+ 



M 2 (M 2 -M 2 -i) 



|a L ||a s |cosO L - (j> s ) 2\a L \ 2 - 3|a s 



11 

T 



± 



|2 2 I |2 2 

a L \ a s - \a s \ a L 



2ollOis&-t 
+ a L a 5 (4|a L | 2 -2|a s | 2 -3) 



(At) 2 



(163) 



(At) 



(164) 



according to the general expression (|58j) . 

Equations (|158p ~ (|164p can be readily generalized to any initial distribution 
of the radiation fields. 



2. Exact solutions 



Let us now proceed to the exact solution of the master equation (|125|) . 
We apply the Laplace transform method. The method is readily applicable 
to nonlinear maste r equations for a variety of nonlinear optical phenomena 
(Refs. [30l [73l f220l | and references therein); in particular, it has been applied 
successfu l ly to diffe rent multiphoton Raman processes in Refs. [73ll75l. [78lll07l 
Hoi . [Til fill. HlH - The solution of (li"25|) for diagonal terms of the density 
matrix p nm (00T) (i.e., for v = fi = 0) was derived by McNeil and Walls [73| 
and then, in a more general form, by Simaan [75| . As the chief result of the 
present work we derive the time-dependence of the complete density matrix 
Pnm{vfJ.T), where the degrees of off-diagonality v, \i are arbitrary. To the best 
of our knowledge, ours is the first derivation of a complete analytical solution 
to the Raman scattering model including depletion of the pump held. 

As usual, we assume that the Stokes and laser beams are mutually indepen- 
dent at the initial time r = tq. Thus, the initial joint distribution p n m{y ^tq) 
is a product of the distributions for the separate beams, 



(165) 



Let us define the coefficient A in terms of the integer- value function [[x]] (the 
maximum integer < x): 



A = 



1 a — v 
- + ! ~^— 



(166) 
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The exact solution of (1125j) under the condition (|165[) . derived in Appendix B, 
reads as follows for A < 0: 



m!(m + jtt)! 



n\(n + v)\ 



1/2 m 



1=0 



(n + l)\(n + l + u)\ 



(m - l)l(m -l + n)\ 



1/2 



£exp[-/(?)Ar] JJ[/(p) - (167) 



whereas for A > it is 







g=0 




is 








m 


(to 




1/2 | 


n 


(n- 


hi/)! _ 





p=0 



(n + l)\(n + I + v)\ 11/2 



(to - l)\(m -l + n)\ 



£exp[-/(g)AT] !][/(?)- /(g)]" 1 

g=0 



p=0 
P^9 



-(l-5 m o) Pn+l( VT 0)Pni-l(P T 0) 



;=a+i 



x 



(n + /)!(n + / + i/)! 



1/2 



(to - /)!(m + 

A l A 



<j=0 <}'=A+1 p=0 

p^g 



p' = A + l 



x(«$ 

+ ( 5 /(<?)/(g') - !) 



/(«)/(?') Arexp[-/(g)Ar] 

exp[-/(g)Ar] - exp[-/(g')Ar] 
/(<?) " /(?') 



where the function f(x) is given by 

f(x) = — [(n + x)(m — x + 1) + (n + x + v)(m — x + /x+ 1], 

The coefficient A can be alternatively defined as 

m — n t n — v 

A 



(168) 
(169) 

(170) 



The solution (|1 6T[) and (|168|) with consequent application of the coefficient A 
of Eq. (I170p reduces, as it should, to the Simaan solution (45) and (49) of 
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Ref. [75| for the diagonal matrix elements [y = fj, = 0); McNeil and Walls have 
also obtained a solution of (|125l) for the diagonal matrix elements (Eqs. (6.5), 
(6.6) and (4.14) in Ref. 73]); however, their solution is not in full agreement 
with Simaan's solution and is not a special case of ours for reasons given by 
Simaan [75| . 

The solution (|167j) and (|168p is very well adapted to numerical analysis; 
nonetheless, it is of a rather complicated form. Our solution of (1125|) can be 
rewritten more compactly. Following the method of Malakyan [116| . we find 
(for details, see Appendix B) 



to! (to + /i)! 



n\{n + v)\ 



1/2 



^Pn+l{ VT 0)P S m-l{^) 



1=0 



(n + l)\(n + I + v)\ 



(to - Z)!(m - I + pS)\ 



1/2 



(171) 



xV exp[-/(g)Ar] [] [/(?) " 



p = 



The differential operator of the d-ih order, X>, is defined as follows: 



T>=(-i) d n 



L Ji df(q r ) 



(172) 



The order d of the differential operator (|172p is equal to the number of pairs 
of mutually equal factors occurring in the product of Eq. (|172[) . f(qi) = f(q[), 
/fe) = /(<?2); •••) f(ld) = fWd)- ^ there are no pairs of equal factors then the 
operator V is defined to be unity (see Appendix B). The solutions (|167[) . (|168p . 
and (|172p represent the chief result of our paper. In Sect. 5.2, in the Raman 
effect model under the parametric approximation we have analyzed, in particu- 
lar, the single-mode solutions for either the Stokes mode or for the anti-Stokes 
mode. For completeness, we give in Appendix C the solution for anti-Stokes 
scattering without the parametric approximation. The degree of off-diagonality 
fi is assumed to be nonnegative (contrary to is); nonetheless, the time depen- 
dence of the complete density matrix p nm (vfiT) is determined by the simple 
relation for the inverse matrix elements: 



(173) 



Thus solutions (|167|) , (|168|) , and / or (|172|) provide an entire specification for all 
measurable properties of the light field under consideration. 

The two- mode (joint) density matrix with elements PnmiypT) (I167P 
and (I168P , enables the calculation of the single-mode (separate) density matrix 
with elements p^ipr) and p^(^r). The Stokes mode matrix elements Pm(t 1T ) 
can be calculated from 



Apt) = p™{vpt) 



(174) 
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and the laser mode matrix elements Pm(pr) can be found analogously with 
the exception that for terms p nm {vpT) with p < the property (|173[) must be 
used. The already mentioned solution of McNeil and Walls [73| corresponds to 
the separate diagonal density matrix (|174p . 

There is yet another manner of expressing the two-mode solutions of the 
master equation (|125|) . p nTO (00r), for any initial distributions, via the density 
matrix elements for the initial number states in the Stokes and laser fields: 



Pnm(OOr) = 



E E 

no— mo— 



P^ mo) {r)p L no {Qr )p s mo (QT ), 



(175) 



where p& ,mo) (00r) is the solution ([157)1 and ([155)1 for p„ m (00r), under the 
initial conditions that the laser field is in the number state |n ) and the 
Stokes mode is in the number state \mo). The weighting functions in (II 75[) 
are arbitrary initial distributions of the laser, p^(Oro), and the Stokes field, 
Pn(® T o)- Here, for brevity, we restrict our considerations to diagonal terms 
(with v = p = 0). Otherwise, instead of pl™' mo) (00T) we would have to use 
the solution pj^m™ ' 1 ' '^ (vpr) and perform two extra summations in Eq. (|175l) 
over v, p. McNeil and Walls [73|] have presented their solution of Eq. (I125[) in 
this manner. Analogously, we can express the single- mode distributions p^(Or) 
(p^(0r)) for arbitrary initial states using the solutions ^ no) (r) (pST o) (t)) for 
the initial photon-number states |no) (|too)). For instance, for the Stokes mode 
solution we apply the formula 



E Pn 

n =0 



no) (r)p^(0r ) 



(176) 



We shall make use of this procedure for the diagonal approximate solu- 
tions (TTM)). 

Having the solutions (I167p . (|168[) . or (|172[) available we can, at least numeri- 
cally, analyze, e.g., the single- and two-mode photocount statistics and quadra- 
ture squeezing. The expectation values in the relations describing squeezing 
and photocount statistics (see Section 4) are readily expressed in terms of the 
density matrix elements Pnmiypr) by way of 



(h k (t)) = n k p n , m (0,0,t), 



(177) 



n.m— 



m k {t)) = ^mVm(0,0,t), 



(178) 



(&l k (t)) = E 



(n + k)\ 



1/2 



Pn,m (M,f), 



(179) 



; w> = E 



(m + fc)! 



TO ! 



1/2 



(0,k,t), 



(180) 
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(a+(t)a+(t)) = VKn + lKm+l^VtM,*), (181) 



(a+(t)a s (t)) = ^[(n + lH'Vll-M). (182) 

n,m 

In Section 3 we defined the s-parametrized quasiprobability distribution 
WW({afc}) and the s-parametrized characteristic function C^({0k}) and ad- 
duced relations between them for any parameter s. Here, we deal with matrix 
elements in Fock basis of the density operator, p nm {vp.T). To achieve consis- 
tency between our analysis of the Raman effect presented in this Section with 
the analysis of Section 5, we shall present some relations between the functions 
({«&}) or C^({/3k}) and the density operator p({ak})- We restrict the 
general formulas for the M-mode fields to our two-mode situations, so that 
({a^}) = {a.L,as,A)- Thes e formulas are in complete analogy with the results 
of Cahill and Glauber (194j for the single-mode case. By virtue of the operator 
T{q.l, ols,a), which is the Fourier transform of the s-parametrized displacement 
operator D^{fl t , 3 ,a) (see Eq. (0), 

f^(a L ,a s ,A) = ^ J D (s \(3 l ,(3 s ,a) 

x exp (a L Pl + a s , A P*s,A ~ cc.) d 2 L d 2 /3 s ,a, (183) 

the density matrix p(cll, as. a) can be obtained from the s-parametrized qua- 
sidistribution W (s) (aL, as, a), (f2T|). 

p{a L ,as, A ) = J yV {s) {a L ,a s ,A) 

xf(- s \a L ,as,A)& 2 a L & 2 as,A- (184) 

The inverse relation, 

W^(a L ,a SA )=TT [p(a L ,as, A )f^(a Ll as,A)] , (185) 

resembles expression (|20|) for the characteristic function C^ s '(0l, 0s. a), which 
is the average value of the displacement operator D^ s \aL, &s,a)- We are inter- 
ested in the relations for the matrix elements of p(oll, as, a)- They immediately 
follow from JTSD and (flgS) : 

p n ,m{v,p) = J yV (S \ a L,US,A) 

x (n,m\f( s \aL,as,A)\n + v,m + p J )d 2 aLd 2 as,A, (186) 



oo oo oo oo 



n—0 m—0 u=—n fj,=—m 

x (n,m\f ( - s ' ) (a L ,as,A)\n + is,m + p), (187) 
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The Fock matrix elements for the two-mode held, 

(n, rn\f lyS \a L ,as.A)\n + v,m + /i) = (n\f^ (a L )\n + v) 

x(m\f<- s \a s , A )\m + p), (188) 

are simply product s of the two single-mode Fock matrix elements given by 
Cahill and Glauber [H3: 

x exp (-J-^KI 2 ) 4*> (*^) WJ , (189) 

where L^\x) is the generalized Laguerre polynomial. The above equations 
show equivalency of the two apparently different formalisms we have been deal- 
ing with: on the one side, the s-parametrized quasiprobability distribution func- 
tions W( s ) obtained within the Fokker-Planck equation formalism presented in 
Section 5, and, on the other side, the density matrix operator formalism under 
discussion in this Section. 

Solutions (|167l) , (|168p , or (|172l) reduce to rather simple expressions in special 
cases, for instance, on the one side, for long periods of time when the laser beam 
is totally depleted, and on the other, for an intense laser beam almost unaffected 
(undepleted) during the process of scattering. We now discuss these two cases. 



3. Long-time solutions 

After a sufficiently long time, the system settles down to a steady state as a 
result of the total depletion of the laser pump. The steady-state solutions can be 
readily deduced from (|167[) and (I168|) . Indeed, in the time limit (r — > oo), the 
nonzero matrix elements p nm (vpoo) must satisfy the condition for the function 
f{x) (|169j) that f(q) = 0, which implies that q — 0. Hence, we have 

Pnmiy^i t = oo) = for n, v ^ Q (190) 

for arbitrary m, m + [i ranging from zero to infinity. All photon-number and 
annihilation operator moments for the laser beam vanish in the time limit 

(n*(oo)) - 0, 

(o£(oo)) =0 for p > 0, (191) 

which reflects the fact that the laser beam is totally depleted. The normaliza- 
tion condition takes the form 

oo 

Pom(OOoo) = I, (192) 

m— 

In the model of hyper-Raman scattering, as was shown by Malakyan jll6l | , there 
intervene in the limit r — > oo the nonzero density matrix elements po m (0fioo), 
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Pim(O^oo), pom(lM°°): an d Plm (~ 1) /-too). Hence, the photon-number moments 
and annihilation operator moments of the laser mode do not vanish, contrary 
to the model of Raman scattering under consideration in view of (|191j) . If we 
assume that initially there are no photons in the Stokes mode, then p^(Oro) = 
/O^(Ooo) (Ref. ^H), which implies that an arbitrary photon- number moment 
(to p (oo)) (with any p) for the Stokes mode in the time limit is identical with 
the corresponding moment for the laser mode, (n p (ro)), at the time r = tq. 



B. Raman scattering without pump depletion 

Compact approximate solutions can be obtained from Eqs. (jl6T[) under the 
condition that the initial laser beam is much more intense than the Stokes 
beam, i.e., (n) >> (m). The depletion of the laser beam and amplification 
of the Stokes beam restrict the validity of this approximation to short evolu- 
tion times t (r << 1). This approximation implies that the density matrix 
elements p nm {v[iT) for A > given by Eq. (|168[) are negligible. Moreover, 
we can simplify the remaining solution (|167|) by setting n ps n ± m. In the 
analysis of the phenomena described by the density /£» nm (^/xr) with small de- 
gree of off-diagonality v (such as quadrature squeezing), we can set n ~ n + v . 
Alternatively, in order to, for instance, investigate phase properties 127Nl29| 



(which require summation over v ranging from zero to infinity) one might as- 
sume that the fluctuations in the laser beam are small in comparison to their 
mean value, i.e., (h) >> y/ ((An) 2 ). Under these approximations the solution 
of (|125[) takes the form 

P nm {v[iT) ps [m!(m + [i)l} 1/2 

m 

J2pn(vTo) P S m -i(»T )l(m - l)l(m - I + pW 1/2 



X 

1=0 

I 

X 



exp[-n(m -q+l + a»/2)At] ]J{q pT 1 - (193) 



q=Q P=o 



Applying the binomial theorem we rewrite (|193j) as 



Pnm{vpr) ps ^ 



m\ I m + fx 

i A i 



1/2 



(e" Ar - 1) 



1=0 

x exp[-n(m + 1 + /i/2)Ar] ^(z/r )p^_ i (^o), (194) 

which, for /i = and v = 0, goes over into Simaan's equation of Ref. [75]. The 
density matrix (|193j) . applied to relations (|177|) - (|182j) . enables the calculation 
of the expectation values and variances for the Stokes mode and the laser mode; 
however, in the latter result of the approximations assumed, we find 

no time dependence of the laser field photon-number moments for a Stokes 
beam initially in a number state containing 



(f[n(r)}) = J2f(n)Pn(OTo)eM-n(m Q + l)}Y,( m * m °)(l-e-n m 



n=0 m— 
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= E/(™v»(° t °) = </( ft )>' ( 195 ) 

The result ([T95k) is valid for any initial number-state Stokes beam, so we 
conclude that the pump beam is time-independent for arbitrary initial Stokes 
beam. The photon-number moments for the Stokes mode calculated from (jl94[) 
are of particularly simple form. For instance, we have 



(m(r)) = (rh) ^exp(nAr)p^(0 



to) 



n=0 



( 



Eexp(nAr)p^(OT )-lJ, (196) 



(m 2 (r)) = ((m 2 ) + 3(m) + 2) ^ exp(2nAr)^(0r ) 

oo 

-3«m) +1)E exp(nAr)p^(Or ) + 1, (197) 



n=0 



(rn(r)n(r)) = ((m) + 1) ^ nexp(nAr)p^(Or ) - (n). (198) 

Equations (|T96| and (fT97|) were obtained by Simaan [75l |. Equation (|196|) is in 
agreement with the Shen relation in Ref. [2h]. Equations (|196|) - (|198p reduce to 
Loudon's results of Ref. [30] for the simpler special case in which no scattered 
photons are excited initially. The sum of the mean photon numbers for the laser 
and Stokes mode (|196p is not a constant of motion, contrary to our former con- 
siderations (|132l) . Nonetheless, in view of the intense laser field approximation, 
the conservation of the total number of photons is at least approximately ful- 
filled. It is easy to find a physical interpretation of Eq. (jl96|) . The first term 
of (|196p describes the amplification of the initial Stokes beam with (rh) pho- 
tons at the time To and can be identified as sensu stricto stimulated Raman 
scattering. The second term of (|196[) corresponds to an amplification of the 
vacuum fluctuations and can be interpreted as spontaneous Raman scattering, 
which occurs even in the case when the Stokes field contains initially no pho- 
tons {(rh) = 0). Note that even in the model of scattering from phonons at 
zero temperature ("quiet" reservoir), spontaneous scattering does take place. 

The coefficients 7s (t) and g\ J s( T )j readily obtained from (|35|) and (|38|) by 
insertion of Eqs. (|196[) - (|198[) . can be explicitly compared to the coefficients 
calculated from other, corresponding relations. Assuming that initially there 
are no photons in the Stokes beam, (rh) — (rh 2 ) — 0, we obtain from small-time 
expansions of the exponential functions in Eqs. (|196[) and (|19T[) the following 

(2) 

simple expressions for the normalized factorial moment 7^ (r): 

7s 2) (r) = 2 7 f + 1 + 2((n 3 ) - (n 2 ) 2 /(n))(n)- 2 Ar, (199) 
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as well as the normalized cross-correlation function g LS (r): 

9 { lh) = 1l ] + ((« 3 ) - (n 2 ) 2 /(n))(n)- 2 Ar/2. (200) 

Equations (I199p and (I200|) can be equivalently obtained from the short-time 
expansions (|134l) and (|139[) . respectively, on omitting the expressions l/(n) and 
(n 2 )/(n) 2 in the terms proportional to At, which are negligible in comparison 
with the terms (n 3 )/(n) 2 and (n 2 ) 2 /(n) 3 . The Simaan approximate relations 

for g^g{r) (32) and 7g (t) (33) in Ref. [7j|, rewritten in our notation (with 
extra —1 in view of ([55)1 and (|3"5|) ). do not reduce exactly to our Eqs. (|199p 
and (|200p . respectively. 

By substituting Eq. (|194[) with f = [i = into (I174p one can obtain solu- 
tion (|176l) . for any initial distribution of the laser mode, with the following 
distribution pm {t) : 

pt o) {r) = exp[-(m + l)n Ar] 

rn 

x£Q(e»°^-l)V ro -,(0r o ), (201) 

1=0 

calculated for the laser field initially in a number state containing no photons. 
In this case the mean ((m(r))) and mean-square number of Stokes photons 
((m 2 (r))), 

(m(r)) = ((m) + 1) exp(n Ar) - 1, (202) 

(m 2 (r)) = ((to 2 ) + 3(m) + 2) exp(2n Ar) 

-3((rh) + 1) exp(n Ar) + 1 (203) 

can be immediately obtained either from (|196l) and (|197[) or from (1201 1) . As- 
suming that the Stokes beam is initially in a coherent state \a), we can perform 
summation in (|201|) which leads to 



./no) 



(r) = cxp[-|a| 2 - n Ar] (l - 
XiFi — m; 1; — |a| 2 (e 



nn At 



-1 



(204) 



where i-Fi is a confluent hypergeometric function. The density matrix elements 
Pm (to) (|204p describe a superposition of coherent and chaotic fields (23l . l22l| . 
This will be more transparent if we rewrite Eq. (|204[) in terms of the average 
number of Stokes photons in the chaotic part, 



(m ch {T)} = exp(n Ar) - 1, 
and the mean number of photons in the coherent part alone, 
(m c (r)) = |a| 2 exp(n Ar). 



(205) 
(206) 
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Then, one obtains, using the Laguerre polynomial L m , the standard form of 
the distribution (PZHf [13, El, [llllll : 

,(no)M _ (^M) m (m c (r)> 



(l + (m cft (r)))i+™ "V l + (m cfc (r)) 

xL ™ f- r nin'+r nff ) ' (207) 

V (m c h(T))(l + (m ch (T))) J 

Similarly, by expressing the relation (|202[) in terms of the mean values (I205P 
and (I206[) it is seen that 

(m(r)) = (m c (r)) + (m ch (r)). (208) 

The general moment of th e p-th order (fh p ) can be found by repeated use of 
the recursion relation |222l | : 

(rrf+Hr)) = (m ch (r))({m ch (r)) + 

+{m c (r))(2{m ch (r)) + + K(r))(m(r)) (209) 

with the help of (|208|) or using the following explicit expression [23|, |42[ : 

(m r (r)) - rlKfc(r)) L r f- /^^, ) ■ (210) 

V {meh{T))J 

The factorial moments of p-th order can readily be calculated from (|209|) 
or (|210p . In particular, the second-order factorial moment reads as follows: 

,(2)/v_ 1 /(™ c (r)) x " 



which takes the minimal value, equal to zero, for the initial time To, since only 
then (rhch) = 0. 

For the Stokes beam initially in a vacuum state |0) the distributions (|204l) 
and (|207[) reduce to the Bose-Einstein distribution 

>°)M = ^(t))"' (212) 

describing a chaotic field (cf. (|T?10 . In this case, in the absence of stimulated 
scattering ((to) = 0), the chaotic field is generated in spontaneous Raman 
scattering as an amplification of the vacuum fluctuations. 

To compare the results for the expectation values of the Stokes mode obtained 
in Sect. 5.2 with the present results, we assume that the laser and Stokes 
beams are initially in a coherent state \ai) and \as), respectively. Performing 
summation in Eqs. (|196[) - (|198[) with the coherent weight function p^(Oro) one 
readily arrives at 

(to(t)) = (|a s | 2 + 1) exp[M 2 ( e A - - 1)] - 1, (213) 
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(m 2 (r)) = (| as | 4 + 4M 2 + 2)exp[M 2 (e 2AT -l)] 

-3(| as | 2 + 1) exp[|a L | 2 (e A ^ - 1)] + 1, (214) 

(m(r)n(r)) = |a L | 2 (|a s | 2 + 1) exp[|a L | 2 ( e AT - 1) + At] - \a L \ 2 . (215) 

Within the Fokker-Planck equation approach under parametric approximation 
(Sect. 5.2) we have obtained Eqs. (|111[) and (|114l) . which can be rewritten, 
using the notation of this Section, i.e. nst = \eL\ 2 jst = \a-L\T, and assuming 
that the mean number of phonons is zero, in the following form: 

(m(r)) = (|a s | 2 + 1) exp(|a L | 2 Ar) - 1, (216) 

(m 2 (r)) = (|a s | 4 + 4|a s | 2 + 2)exp(2|a L | 2 AT) 

-3(|a s | 2 + l)exp(|a L | 2 Ar) + l. (217) 

We also note that 

(m(T)n(T)) = (m(r))(n). (218) 

For short times of evolution, At << 1, and intense pump beams, |ai| 2 >> 1, 
Eqs. (|2T3|) . (|2T4]) . and (|2l"5j) go over into Eqs. (|2T6|) . (prTj) . and (j2l8|) . respec- 
tively. Indeed, the short-time expansions of (|213p and (|214[) are: 

(m(r)) = |a s | 2 + |a L | 2 (l + |a s | 2 )Ar 

+ |a i | 2 (|a i | 2 + l)(l + |a s | 2 )^, (219) 

(m 2 (r)) = |as| 2 (l + |as| 2 ) + |a L | 2 (l + 5|as| 2 +2|a5| 4 )Ar 

+ |a L | 2 (|a L | 2 + 1)(5 + 13|a s | 2 +4|a s | 4 )^^, (220) 

whereas Eqs. (|216[) and (|217p obtained within the formalism of Sect. 5.2 reduce 
to 



(m(r)) = |« s | 2 + |a L | 2 (l + |a s | 2 )Ar 

+I«l| 4 (i + |«5| 2 )^?-, (221) 



(m 2 (r)) = |a s | 2 (l + M 2 ) + M 2 (l + 5|a s | 2 + 2|a s | 4 )Ar 

+ |a L | 4 (5 + 13|a s | 2 + 4|a s | 4 )^^, (222) 

respectively. It is seen that for high intensity of the pump field, (|219j) goes over 
into (|2"2Tj) . and ([22TJ1) into (1222]) by setting \a L \ 2 (\a L \ 2 + 1) w |a L | 4 . Hence, the 

factorial moment "t^\r), 

7 f(r) = 2|a L | 2 |a s r 2 Ar- [|a L | 2 (3 + |a s | 2 ) 

-|a s | 4 - 5|a s | 2 - 2l |« L | 2 |a s r 4 (Ar) 2 , (223) 
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calculated from (|219l) and (|220|) in the case of nonzero as and an intense pump 
beam, goes over into 

7 f (r) = 2|a L | 2 |a s |- 2 Ar - (3 + \a s \ 2 )\a L \ 4 \a s \-HAr) 2 (224) 
obtained from (|22ip and (|222l) . If the initial held contains no Stokes photons, 

(2) 

we obtain the following factorial moments 7^ (t): 

7 f (t) = 1 + 2\a L \- 2 + 2Ar + (2 + ^\a L \ 2 ^j (At) 2 , (225) 

7 f(r)=l, (226) 
within the formalisms of this Section and Sect. 5.2, respectively. The diffcr- 

(2) 

ences between the factorial moments j s (r) are more pronounced in the case 

as = since the expansion of (m(r)) and (m 2 (r)) correct to the third order 

in r is required in the derivation of (|225l) . The interbeam degree of coherence 
(2) _____ 

9ls( t ) PHI) , as expected, vanishes for the model of Sect. 5.2. The short-time 

(2) _____________ _____ 

expansion of g LS ( T ) obtained from Eqs. (|213[) - (|215[) . for as ^ 0, is 

g (2) s (r) = (l + |a s |- 2 )Ar 

-(1 + |as| 2 )(2|a L | 2 - |as| 2 )|a5r 4 ^, (227) 
otherwise, as — 0, we get 

9ls(t) = Wl\- 2 + \at + i«n) + 3)^. (228) 

It is seen that the approaches of Sects. 5.2 and 6.2 give similar predictions for 
the Stokes beam. 

The evolution of the photon-number moments is demonstrated in Figs. 2 
and 3: (m(r)) calculated with Eq. (|22ip is depicted by solid line C or with 
(f2T9|) by solid line D in Fig. 2; (m 2 (r)) obtained from Eq. ((222)) is given by 
solid line C or from (|220j) by solid line D in Fig. 3. No time dependence of 
(n(r)) and (n 2 (r)) is observed for the results of this Section and Sect. 5.2, i.e., 
we obtain straight lines C and D in Figs. 2 and 3. Similar notation is used in 

(2) (2) 

Figs. 4-6 for the normalized factorial moments j s (t) (Fig. 4), Jl (t) (Fig. 

(2) 

5), and the degree of interbeam coherence g^s( T ) (Fig- 6). We have chosen 
rather small initial numbers of laser photons (|«l| 2 = 2) for numerical reasons. 
In this case, the factorial moments, calculated from (I223p . (|225p . (|227p . and 
(|228p . differ significantly from the exact numerical results. So we omit them 
(curves D) in Figs 4 and 6. 

VII. CONCLUSIONS 

In this paper we have presented in a systematic way quantum-statistical the- 
ory of the standing- wave Raman scattering within the s-parametrized quasidis- 
tribution formalism in Section 5 and the density-matrix formalism in Section 
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6. Particular attention has been paid to quantum properties of light such as 
squeezing and sub-Poissonian photon-counting statistics. 

In Figs. 2-9 we compared various statistical moments obtained (i) from 
numerical calculations utilizing the exact solutions without the parametric ap- 
proximation - (Eqs. (|167p . (I168p . and/or (|172[) derived in Appendix B), (ii) 
from the short-time solutions of Sect. 6.1.1, (iii) from the solutions obtained 
in Sect. 5.2 within the framework of the FPE approach under the parametric 
approximation, and (iv) from the approximate solutions derived in Sect. 6.2 
within the density-matrix formalism. 

In Figs. 2,3,7 and 8 we demonstrated that the initial and short-time behav- 
ior of the approximate functions (curves B,C, and D) is consistent with the 
exact evolution (curves A). We have shown analytically that our expressions 
for the Stokes scattering formalisms presented in Sects. 5.2, 6.1.1, 6.1.2, 6.2 
are equivalent for short times and high initial intensities of the pump field. 
Nevertheless, it is seen that the equations derived in Sect. 6.1.1 give the best, 
whereas those derived in Sect. 6.2 give the worst approximation to the exact 
solution of Sect. 6.1.2 for small initial intensities of the laser field. 

We showed in Fig. 4 (curve A) that the Stokes mode photon-number fluctua- 
tions vary from initially chaotic to Poissonian in asymptotics (for (m) = and 
(h) = |ai| 2 ) or from Poissonian, through super-Poissonian, to Poissonian for 
large times (if (m) = |«s| 2 7^ and (n) — |«l| 2 )- The asymptotic behavior of 
is ( T ) * s consistent with our predictions in Sect. 6.1.3. The short time behav- 
ior of 7c; 2 ■*("!") for hyper-Raman scattering jll6l | is similar to that presented in 
Fig. 4 for Raman scattering. However, for long times the Stokes hyper-Raman 
photon- number fluctuations become sub-Poissonian for reasons given in Sect. 
6.1.3. 

In Fig. 5 we demonstrated that the normalized factorial moment for the 
(2) 

laser mode, j L (t), varies from initially Poissonian to super-Poissonian. The 

differences in 7^ '(t) between Fig. 5a (spontaneous Stokes scattering) and Fig. 
5b (stimulated Stokes scattering) are only quantitative. We note that for hyper- 
Raman scattering the photon- num ber fluctuations in the initially coherent laser 
mode become sub-Poissonian [ll6| . The time behavior of the interbeam degree 
of coherence was presented in Figs. 6a and b. Curve A in Fig. 6a coincides 
with the Simaan exact solution [75j . Sub-Poissonian statistics in the compound 
laser-Stokes mode is observed. 

It is thought (see, for instance Ref. (22|, p. 192) that the Simaan approach [75| 
is the most rigorous application of master equations to the Raman problem. 
However, Simaan's solution is restricted to the diagonal matrix elements in 
number representation. Only these terms are needed to obtain the mean pho- 
ton numbers and their higher moments, whereby the photon correlation effects 
can be investigated. To investigate squeezing properties and phase correla- 
tions it is necessary to obtain the off-diagonal elements of the density matrix. 
We generalized the solution of the master equation obtained by McNeil and 
Walls [73| and Simaan [75| to comprise all the off-diagonal matrix elements as 
well. Our derivation of the complete density matrix represents the main result 
of this paper. 

Our intention was to cite an extensive literature related to our Raman scat- 
tering approaches. Nevertheless, we realize that the cited literature is not com- 
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plete. We include only those references that are most relevant for the purposes 
of our article. 
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Appendix A: Solution of generalized Fokker-Planck equation 



Here, we give a simple formal solution of the generalized FPE (fBTj) for the 
quasidistribution 1 '(ul, as, <%A, t) (Q-functions) as well the correspond- 
ing characteristic function C^ 1 '{Pl, As, Pa, t) - the solution of the simplified 
equation of motion (|63p. We choose antinormal order (s — —1) to avoid the 
problems of existence of the quasiprobability distributions and to reduce the 
FPE (fodj) containing the third-order derivatives (for s ^ ±1) to a second-order 
FPE. It is by no means easy to find an exact solution of (|5Tj) or even for 
particular orders, because the drift coefficients are not linear and the diffusion 
coefficients are not constant. An often employed method to solve problems of 
this kind is to assume that the fluctuations of the radiation fields are small 
compared to their mean v alues ; i.e., the q uasidistribution describing the fields 
is sharply peaked [H H [l87l [HI Hi- This will be the case for suitably 
chosen input state and the initial output states. Under these restrictions we 
can rewrite our FPE (1611) related to antinormal order in the linearized form 



d d 
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where the coefficients Dm, Dm, Ck, £k for k,l = L,S,A are defined by (jSlIj) at 
the i nitia l moment to. Equation (IA1|) is the generalization of the FPE given in 
Ref. |217l | for the case of nonzero 7s and jas- It is seen that the Raman effect 
under t he ap proximations applied can be treated as an Ornstein-Uhlenbeck 
process jl88l |. The FPE (|A1|) can be solved exactly by various techniques; 
see, e.g., Ref. [H3. For instance, using the inverse Fourier transform (|22[) . 
one can transform the FPE (|A1[) into the corresponding equation of motion for 
the characteristic function C^ _1 ^(/3l, Ps, Pa, t), which is a first-order differential 
equation. The method of characteristics applied to the latter equation leads to 
the solution 



exp 



B 



(-1) 



k=L,S,A 



> r [-Ci(t)(3l+c.c.)+{Ml{t)-cc.) 



+[D L s(t)Pl0* 3 + D L sPL,P* s 

+ D LA {t)PlP*A + D L APLP*A 

+D SA (tW A + DsaPsPa + cc] } 



(A2) 



The angle brackets mean averaging over the complex amplitudes (fe = 
L, S,A) with the initial distribution W^ 1 ^ (aL, as, a a, to)- They represent 
the influence of the initial photon statistics of the pump and scattered fields on 
the evolution of the system. The solution of Eq. (IA1|) can be readily obtained 
by applying the Fourier transform (|2"Tj) to solution (| A2|) . and has the form of 
a shifted seven-dimensional (including time) Gaussian distribution involving 
correlation between the radiation fields, 

W^\ ai ,a 2 ,a 3 ,t) = /-^expl-^- 1 ))- 2 ^ 



3=1 



3=1 k=j+i 
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where, for simplicity, we have identified the subscripts in a\ — ul, oli = as, 
a3 = &A, 6 = £l, 6 = £s, 6 = £a- The functions E{ 1 >, ifo , and 
which are time-dependent, are connected with the functions 1 \ Cfc, Dki, 
Dki appearing in (| A2|) in a manner similar to (|74|) and (|75|) . respectively. We 
do not adduce explicit formulas for the coefficients listed, since solutions (|A2j> 
and (|A3|) serve only as an example of how one can deal with Eqs. (J6TJ) and (|63|) . 
The validity of solutions (|A2[) and (|A3|) is restricted by strong approximations, 
which are actually equivalent to the parametric approximation and the short- 
time approximation. 



Appendix B: Solution of master equation in Fock basis 

Here, we solve the equation of motion (I125[) . To eliminate the square root 
appearing in Eq. (|125[) for off-diagonal terms, it is convenient to introduce the 
transformation 



n\(n + v)\ 11/2 



to! (to + /i)! 



>mt), (Bl) 



where the degree of off-diagonality [i is restricted to nonnegative integers, 
whereas the degree v is > —to. On insertion of (|B1[) into (|125p . the equa- 
tion of motion for the transformed matrix elements ip n m{v^,T) takes the form 

-0nm(^ r ) = -\[ n { m + l ) + (n + v)(m + H + l)]tf) nm (vfJ,T) 

+ 1pn+l,m-l(v^T). (B2) 



We apply the Laplace transform method to (|B2[) . which readily leads to the 
solution 

m I 

V'nmC^Ms) ^^^n+l.m-l^flT^Ylis + .fip)}^ 1 (B3) 
1=0 p=0 

for ?/; nm (^/iT), the Laplace transform of ip n m(v '^ T ) ■ The function f(p) occurring 
in (|B3I) is given by (|169l) . If there are no equal terms among the elements of 
the set /(0), /(l), f(l) the inverse transform, after retaining the /0„ m (^/ir) 
notation, yields (|167l) . If there are repeated elements in the denominator of 
(|B3|) . the inverse transforms will involve convolutions. We apply two ge nera l 
procedures essentially equivalent to that of Simaan [75[ and Malakyan [116| . 
It is convenient to split (|B3I) into two terms as follows: 

a i 

1=0 p=0 



n-\-l,m—l 



l=\+l 
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A I 

xflta + ZCp)]- 1 II [ a + /(p')] _1 (B4) 

p=0 p'=A+l 

with A defined by (|166p (or equivalently by (j 1 70[> ) . Let us note that a parabola 
/(g)=const takes its maximum value for go = (2m- — 2n + \i — v + 2)/4. This 
value, qo or better A, the maximum integer < q$, can serve as a criterion to 
split (IB3|) in such a manner that a convolution theorem can be easily applied. 
The first term in (|B4j) has no mutually equal factors in the denominator, so the 
inverse Laplace transform has the form of (|167l) with the proper upper limit 
of summation. The denominator of the second term of (IB4[) contains repeated 
factors, which are separated, so that we can readily apply the convolution 
theorem finally obtaining the solution (|168j) . Equations (|167j) and (|168p have 
a rather complicated structure. We can rewrite (|167[) and (I168[) in a more 
compact form. If we assume that there is only one pair of equal factors among 
the elements of the set /(0), /(l), fit), i.e., if 

V = M) A * 

7 = 0,-..,! 

9i,«ie{o,...,i} «' 4 "-«i 
then we can express the solution (|B3|) for ip nm (ufir) as 

m Z 
^ml^s) = ^Vn+i,m-i( i/ / x ' r o)[s + /(Q'l)]" 2 JJ [s + /(p)] -1 . (B5) 

i=0 p=° 



The inverse transform of (|B5|) is 

m 

n+Z,m— Z (z^r ) 

1=0 

{/ z 
£ exph/^Arint/W-Zte)]- 1 



q=0 p=0 



( 



Ar- ^ [/(fc) - fiqiT 1 ) exp[-f( qi )AT] 

k=0 
p=0 

which is a derivative of the solution (|167[) over /(gi) with extra minus, 

m 

( 



(=0 

/ 



a/(gi) 



(B6) 



£ exp[-/(g)Ar] J] [/(p) " ( B7 ) 
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In the case of d equal pairs, i.e., f(qi) = f(q[),—,f(<ld) = f{q'd), the Laplace 
transform solution (IB3I) can be rewritten as 



1=0 



n^+z^')]- 2 n i s +f(p)} 



p=Q 



finally leading to 



1=0 



(B8) 



\ df(q r ) J 

]T exp[-/( 9 )Ar] n [/(P)-/(?)] _1 . ( B9 ) 



or equivalently to the solution (| 1 72[) with the d-th order differential operator 

v (HZU). 



Appendix C: Solution for anti-Stokes scattering 



In Section 6 we have given an analysis of Stokes scattering. For complete- 
ness, in this Appendix, we present the solution describing the anti-Stokes effect 
including laser depletion, but neglecting the Stokes generation and assuming 
that the reservoir is 'quiet', i.e., 



IS = JAS = ISA 

(n v ) = 0. 



0. 



(CI) 



Under these conditions the master equation (|18[) in Fock representation is 
d 1 

— p n7n (vfiT) = --[(n + l)m+(n + v+l)(m + y)]p nm (vfiT) 
or 2 

+ [n(n + v)(m + l)(m + (i+ 1)] 1/2 

Pn— l,m.+l 

far), (C2) 

where, for brevity, we have set — n, n' L = n + v (as in Section 6) and 
riA — m, n' A = m + fi. If we define A as follows: 



A 



n — m + 1 v — (i 



2 4 
then the solution of (IC2I) for A < becomes 



(C3) 



n\(n + v)\ 



to! (to + 



1/2 n 

^Pn-liv^pt+li^o) 
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(m + l)l(m + l + fi)\ 



1/2 



(n-l)\(n-l + v)\ 
i i 
x ^exp[- 5 (g)Ar] \[[g{p) - g(q)}-\ 

g=0 



(C4) 



p=0 



whereas for A > it becomes 



nl(n H 


-i/)! " 


1/2 | 


to!(to - 


f-/i)!_ 





(m + /)!(m + / + /i)!" |1/2 



{n-l)\(n-l + v)\ 
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i 



x^ex P [- 9 (g)Ar]rj[ 5 (p)- 9 (g)] 1 
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x 
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A i A 
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q=0 q'=\+l P=° 

P^Q 



p' = A + l 

p'/g' 



x ^s(g)9(g') Arex P[ _ S , (Q')AT] 

exp[- 5 (g)Ar] - exp[-g(q')Ar] 



g(i)g{i') 



1) : 



5(9) - g(q') 



1, 



5(2:) = — [(m + cc)(n — x + 1) + (m + x + /i)(n — at + f + 1]. 



with 



Alternatively, we can express solution (|C4I) and (|C5|) as 
n\(n + v)\ ~\ 1 ■ 



(C5) 



(C6) 



m!(m + /i)! 



(=0 



(m + l)l(m + l + 



(n — l)\(n — l + v)\ 



1/2 



xf) ^2 exp[- ff ( g )Ar] J] fo(p) - 



(C7) 



using the differentiation operator V given by (|172l) . 
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